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INTRODUCTION 


The  puipose  of  this  technical  report  is  to  document  experimental  and  modeling  results 
obtained  during  FY  89  and  FY90  on  the  High  Pressure  Combustion  Kinetics  (E3)  and 
Combustion  Sciences  (Q2)  projects.  This  work  was  performed  in  the  Combustion  Research 
Laboratory  at  the  Astronautics  Laboratory  (AFSC).  In  general,  the  experiments  involved  the 
spectroscopic  study  of  propellant  flames.  The  modeling  work  involved  chemical  kinetic 
modeling  of  solid  propellant  flame  chemistry.  Previous  work  can  be  found  in  References  1-8. 

The  next  section  ("Laser-Induced  Fluorescence  Measurements")  reports  attempts  to 
apply  laser-induced  fluorescence  (LIF)  measurements  of  various  molecules  to  solid  propellant 
flames.  "Hydrazine  Diffusion  Flame  Experiments"  presents  initial  results  obtained  in  the  study 
of  hydrazine  diffusion  flames  at  low  pressure.  The  last  section  ("Chemical  Kinetic  Modeling") 
presents  the  chemical  kinetic  modeling  results.  The  buildup  of  the  hydrazine  burner  occurred 
in  projects  G5  and  Q3.  Future  results  obtained  from  this  burner  will  be  reported  under  the  UV 
Mechanisms  project  (N3). 

The  composition  of  the  propellants  discussed  in  this  report  is  listed  in  Table  1.  Further 
information  can  be  found  in  Reference  5. 


Table  1.  Propellant  Formulations. 


HMX1 

HMX2 

API 

AP2 

AN1 

DB1 

DB2 

Ingredient,  wt  %  (approx) 
HMX  (200/20  pm) 

73 

80 

0 

0 

0 

0 

0 

polyester  binder 

10 

20 

0 

0 

0 

29 

27 

TMETN 

17 

0 

0 

0 

12 

0 

0 

NG 

0 

0 

0 

0 

0 

71 

67 

AP  (400/20/2  pm) 

0 

0 

87 

68 

0 

0 

6 

Polybutadiene  binder 

0 

0 

13 

16 

0 

0 

0 

AN 

0 

0 

0 

0 

16 

0 

0 

GAP  binder 

0 

0 

0 

0 

16 

0 

0 

Aluminum  powder 

0 

0 

0 

16 

0 

0 

0 

HMX=cyclotetramethylene  tetranitramine  C^gNgOg,  TMETN=trimethylolethanetrinitrate 
C5H9O9N3,  NG=nitroglycerin  C3H5N3O9,  AP=ammonium  perchlorate  NH4CIO4, 
AN=ammonium  nitrate  NH4NO3,  GAP=glycidyl  azide  polymer.  The  "polyester"  binder  in  the 
DB  propellants  is  a  polyethylene  glycol-based  polymer  ("PEG"),  while  the  HMX  binder  is 
based  on  polydiethylene  glycol  adipate  ("PGA"). 
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LASER-INDUCED  FLUORESCENCE  MEASUREMENTS 


OH  UF  TEMPERATURE  MEASUREMENTS 

In  recent  OH  LIF  measurements  at  the  Astronautics  Laboratory,  a  detection  system 
employing  a  much  smaller  bandpass  coupled  with  a  double  spectrometer  (Spex  1404)  has  been 
used  to  increase  the  light  rejection  capabilities  of  the  experiment,  compared  with  earlier  work 
[6,7].  The  bandpass  was  decreased  from  20  nm  to  0.45  nm,  as  shown  in  Figure  1.  Wide 
detection  bandpass  is  desirable  because  it  greatly  simplifies  the  calculation  of  a  rotational 
temperature  [9].  This  decrease  in  bandpass  increases  the  light  rejection  capability  of  the 
detection  system  at  the  expense  of  ease  of  calculating  a  rotational  temperature  from  the  data. 
The  narrow  bandpass  was  centered  on  the  OH  1,1  Ri  bandhead  so  that  LIF  from  the  OH  1,1 
Rl  1-Ri  10  transitions  was  collected.  It  was  hoped  that  this  detection  bandpass  would  increase 
light  rejection  while  providing  an  accurate  sample  of  the  LIF.  Successful  OH  LIF 
measurements  were  made  using  this  system,  as  shown  in  Figures  2  and  3.  In  Figure  2,  an  OH 
LIF  profile  in  an  atmospheric  pressure  (N2  atmosphere)  flame  of  an  87%  ammonium 
perchlorate,  13%  polybutadiene  binder  propellant  is  shown.  The  profile  was  measured  by 
allowing  the  propellant  to  bum  down  through  the  laser  beam/collection  volume  [4,5],  The  ' 
transmission  of  the  laser  across  the  surface  is  also  shown,  giving  a  relatively  accurate  location 
of  the  surface.  The  "spikes"  in  the  LIF  data  below  the  surface  come  from  reflections  of  the 
laser  off  the  side  of  the  propellant  strand.  The  corresponding  experiment  with  the  laser  detuned 
from  an  OH  resonance  is  shown  in  Figure  3,  demonstrating  the  capability  of  the  detection 
system  to  differentiate  between  OH  LIF  and  scattering  at  the  laser  wavelength.  This 
discrimination  didn't  occur  in  the  "wide"  bandpass  experiments  [6,7],  The  experiment 
corresponding  to  Figure  2  was  attempted  in  the  HMX1  propellant  at  18  atm,  but  no  detectable 
OH  LIF  signal  was  measured.  This  is  perhaps  not  surprising  since  the  equilibrium  OH 
concentration  in  the  AP  propellant  flame  is  15  times  greater  than  that  in  the  HMX1  propellant 
flame  (0.014  vs  0.0009  mole  fraction),  indicative  of  the  dramatically  different  stoichiometries 
in  the  two  propellant  flames  (calculations  made  with  NASA  Lewis  thermochemical  code). 
Quenching  is  probably  much  more  important  at  the  higher  pressure  also,  decreasing  the  LEF 
signal  strength.  The  "good  news"  is  that  no  scattering  was  seen,  indicating  that  the 
measurement  can  be  successful  in  a  high  pressure,  highly  scattering  environment,  provided  that 
a  sufficient  amount  of  OH  is  present  For  comparison,  the  equilibrium  OH  mole  fraction  in  a 
pure  HMX  flame  at  18  atm  is  0.017,  indicating  that  such  an  OH  LIF  temperature  measurement 
in  a  strand  of  pure  HMX  might  be  successful. 

If  one  has  two  profiles  such  as  the  one  shown  in  Figure  2,  with  the  only  difference 
between  the  two  being  the  OH  transition  excited,  then  it  is  relatively  simple  to  calculate  a 
rotational  temperature  [6,10].  This  was  done  in  the  AP  propellant  flame  using  the  OH  1,0  Ri4 
and  Ri  10  transitions.  Note  that  the  directly  coupled  OH  1,1  Rj4  and  Ri  10  lines  (312.395  and 
312.247  nm,  respectively)  are  within  the  detection  bandpass.  Then,  with  the  two  excitation 
lines  denoted  by  subscripts  1  and  2:  kT/hc=(Erl-Er2)/ln[(In/g"B)i/(Ifi/g"B)2]  in  the  unsaturated 
case  and  kT/hc=(Eri-Er2)/ln[[Ifl(l/g"+l/g')]i/  [Ifl(l/g"+l/g')]2}  for  the  saturated  case.  In 
these  equations.  In  is  the  collected  LIF  intensity,  Er  is  the  rotational  energy  of  the  transition,  B 
is  the  absorption  coefficient  and  g  is  the  degeneracy  [6].  The  reproducibility  of  the  data  was 
only  fair,  so  three  experiments  with  each  excited  transition  were  averaged  and  the  resulting  OH 
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rotational  temperature  profile  is  shown  in  Figure  4  for  the  unsaturated  case.  For  highest 
accuracy,  the  data  should  be  corrected  for  partial  saturation  [6,1 1].  In  the  partially  saturated 
regime,  the  rotational  temperatures  yielded  by  the  main  and  satellite  branches  differ  when  either 
total  saturation  or  unsaturation  is  assumed.  Since  this  correction  was  not  made,  the 
temperature  profile  in  Figure  4  is  only  a  rough  estimate  at  this  point.  To  calculate  a  more 
accurate  rotational  temperature  from  the  narrow  bandpass  data  will  require  background 
experiments  in  a  flame,  ideally  a  flame  with  a  similar  pressure  and  stoichiometry  to  the 
propellant  flame.  In  some  preliminary  work  in  a  CH4/N2O  flame,  an  excitation  scan  was 
performed  with  the  narrow  detection  bandpass.  The  result  is  shown  in  Figure  5.  Rather  than 
the  straight  line  found  with  wide  detection  bandpass  [6,10],  the  Boltzmann  plot  shows  some 
substantial  scatter.  The  temperature  is  also  higher  than  that  found  with  wide  detection 
bandpass.  Thus,  accurate  temperature  measurements  will  require  calibration  in  a  flame  with  a 
known  temperature  profile.  In  summary,  then,  it  has  been  demonstrated  that  OH  LIF 
temperature  measurements  can  be  made  with  this  type  of  detection  system,  although  the  OH 
concentration  must  be  relatively  high.  OH  absorption  temperature  measurements  have  also 
been  reported  in  high  pressure  solid  propellant  flames  [12],  as  have  CARS  temperature 
measurements  [13].  Both  of  these  techniques  have  lower  spatial  resolution  than  the  potential 
possible  for  OH  LIF  temperature  measurements,  but  have  yielded  more  successful  results  thus 
far. 


Figure  1.  Bandpass  comparisons  for  OH  LIF  temperature  measurements.  The 
narrow  bandpass  is  312.25±0.22  nm,  encompassing  the  OH  1,1  Ri4-Rill  lines. 
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Figure  2.  Results  of  OH  LIF  "bumdown"  experiment.  AP  propellant,  1  atm  N2, 
0.2  mJ/pulse  in  flame,  50  Hz,  OH  1,0  RilO  excitation  (281.537  nm),  slit:  200  pm  X  1  cm 
projected  1:1  on  flame  (-spatial  resolution),  single  pulse  (no  averaging),  1.7  mm/s  propellant 
bum  rate.  Data  offset  from  zero  on  Y  axes  for  clarity. 
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Figure  3.  Same  as  Figure  2,  except  laser  detuned  from  OH  resonance  lines 
(281.11  nm). 
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Figure  4.  AP  propellant  temperature  profile.  Obtained  for  unsaturated  case  with  Ri4 
and  RilO  profiles  (e.gM  Figure  2).  The  average  of  three  experiments  for  each  excitation  line 
was  used.  The  adiabatic  flame  temperature  for  this  propellant  is  about  2900  K. 


Figure  5.  Results  of  excitation  scan  (Boltzmann  plot)  in  CH4/N2O  flame  using 
"narrow”  detection  bandpass.  0.4  mJ/pulse  in  flame,  laser  scanned  from  281.1  to  281.6  nm 
(OH  l,0Ril-Ril0). 
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CN  UF  MEASUREMENTS 


Several  CN  LIF  excitation  spectra  have  been  reduced  obtain  a  rotational  temperature. 
This  may  be  of  interest  as  an  alternative  temperature  measurement  in  the  reaction  zone  of  flames 
although  OH  would  be  better  for  general  use  since  OH  is  found  throughout  flames,  while  CN 
is  found  only  in  the  "reaction  zone."  Relatively  extensive  CN  LIF  experiments  in  solid 
propellant  flames  have  been  performed  at  AL  [4,5]  although  the  signal  levels  would  be 
marginal  at  best  for  temperature  measurements  similar  to  those  discussed  below.  In  any  case, 
both  the  propellant  and  excitation  scan  measurements  were  made  by  exciting  the  CN  0,1 
transition  at  ~421  nm  and  detecting  the  CN  0,0  band  at  —388  nm,  as  shown  in  Figures  6-8.  In 
Figure  7,  the  CN  0,0  LIF  spectrum  from  excitation  of  the  CN  0,1  P  branch  bandhead  at  421.6 
nm  is  shown.  In  Figure  8,  the  CN  0,1  excitation  scan  is  shown,  with  the  detection  bandpass 
of  0.09  nm  approximately  centered  on  the  388.3  nm  CN  LIF  P  branch  peak  shown  in  Figure 
7.  These  measurements  were  performed  in  a  conical  premixed  CH4/N2O  flame  above  a  0.2 
inch  tube  [5],  at  the  tip  of  the  "inner  cone",  where  the  CN  LIF  signal  is  maximum.  The  result 
of  a  brief  examination  of  the  CN  LIF  intensity  vs  laser  power  is  shown  in  Figure  9.  Thus,  the 
excitation  scan  of  Figure  8  is  apparently  partially  saturated.  A  problem  with  this  scheme  is  the 
possibility  of  self-absorption  of  the  emitted  radiation  in  the  0,0  band.  Significant  self¬ 
absorption  has  been  seen  in  a  50  torr  CH4/N2O  premixed  flat  flame  in  this  laboratory  using  0,1 
excitation/0,0  detection. 

To  obtain  a  temperature  from  Figure  8,  first  the  CN  lines  must  be  identified.  Only 
limited  literature  data  was  found  for  the  0,1  band  (N=l-29,  doublets  not  resolved  [15]),  so  the 
CN  0,1  lines  were  calculated.  This  was  done  by  creating  a  "synthetic"  CN  0,1  spectrum  using 
standard  techniques  and  data  [14-17].  The  program  is  a  modification  of  the  program  listed  in 
Reference  5,  with  the  data  listed  in  Tables  2  and  3.  The  model  data  matched  literature  CN  line 
data  fairly  well,  as  is  shown  in  Figures  10  and  11.  The  model  and  experimental  CN  LIF  data 
were  then  compared  and  line  identification  was  relatively  easy,  as  shown  in  Figures  12, 13  and 
14.  The  peak  heights  were  then  (roughly)  measured,  and  Boltzmann  plots  were  constructed  as 
shown  in  Figures  15  and  16.  The  same  process  was  performed  for  an  excitation  scan 
performed  at  lower  power,  with  the  Boltzmann  plots  shown  in  Figures  17  and  18.  The 
resulting  temperature  of  about  1700  K  (only  the  "unsaturated"  temperatures  were  reasonable) 
were  somewhat  lower  than  the  2600  K  measured  in  the  same  flame  with  OH  LIF  [6],  which 
seems  possible  since  the  CN  excitation  scans  were  performed  in  the  region  of  the  flame  with 
maximum  CN  signal  (the  "reaction"  zone),  while  the  OH  LIF  excitation  scans  were  performed 
in  the  region  above  the  reaction  zone.  However,  modeling  and  experimental  measurements  of 
CN  in  flames  in  References  7  and  11  and  in  the  modeling  section  of  this  report  show  the  CN 
peak  occurs  where  the  flame  temperature  first  reaches  its  maximum.  Thus,  it  is  also  possible 
that  the  excitation  scans  are  yielding  temperatures  that  are  substantially  lower  than  the  correct 
value.  Perhaps  the  fact  that  the  laser  is  exciting  v"=l  is  contributing  to  this  error.  It  would  be 
interesting  to  determine  the  CN  rotational  temperature  in  a  flame  with  a  known  T  distribution. 
Since  the  LIF  is  collected  from  the  0,0  band,  there  is  a  possibility  that  self-absorption  is  also 
important  The  partial  resolution  of  the  two  spin  components  also  leads  to  errors  in  the  data. 
The  partial  resolution  of  the  CN  lines  was  corrected  for  when  reducing  the  excitation  scan  data 
by  calculating  a  correction  factor  based  upon  how  well  resolved  the  peaks  were.  For  totally 
unresolved  peaks  (N"£20),  the  factor  was  1.  Totally  resolved  peaks  had  a  factor  of  0,  with  the 
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partially  resolved  peaks  in  between.  The  peak  height  was  then  divided  by  1+comcction  factor. 
For  example,  the  height  of  the  Pi  10  +  P2  10  unresolved  doublet  was  divided  by  2  to  obtain  the 
heights  of  the  individual  Pi  10  and  P2  10  lines,  while  the  Pi  37  peaks  were  divided  by  1.3  to 
account  for  the  contribution  of  one  peak  to  the  height  of  the  other.  Another  obvious  conclusion 
is  that  the  R  branch  would  be  better  for  temperature  measurements,  since  there  is  no  bandhead 
to  confuse  the  measurement  Note  that  the  structure  of  the  CN  0,1  P  branch  is  inconvenient  for 
temperature  measurements  for  another  reason:  the  rotational  levels  with  the  highest  population 
(Fig.  19)  coincide  with  the  bandhead.  As  can  be  seen  by  the  laser  power  vs  wavelength  curve 
in  Figure  8,  however,  the  dye  used  in  these  experiments  (Coumarin  120)  would  not  be  useful 
for  such  a  measurement.  Bis-MSB  would  be  better  [5].  In  reducing  the  data,  the  line  strength 
was  calculated  as  in  Herzberg,  rather  than  using  tabular  Einstein  coefficients  as  in  OH 
temperature  measurements  [17]. 

Another  experiment  was  performed  that  included  more  CN  lines.  The  new  lines  are 
illustrated  in  Figures  20  and  21.  The  lines  were  easy  to  assign,  as  shown  in  Figure  22. 
However,  the  Boltzmann  plot  still  yielded  lower  temperatures,  as  shown  in  Figure  23. 
Another  possible  source  of  error  is  the  detection  bandpass  used  (0.09  nm  for  Figures  15,16 
and  0.22  nm  for  Figures  17,18).  Assuming  that  the  bandpass  is  centered  to  include  the 
maximum  number  of  CN  lines,  a  bandpass  of  0.09  nm  would  include  N=20  to  36  P  branch 
lines,  while  a  bandpass  of  0.22  nm  would  include  N=15  to  41.  Thus,  the  bandpasses  for 
Figures  15-18  includes  the  directly  excited  lines  at  high  rotational  energy,  but  not  at  low 
rotational  energy.  Thus,  the  two  sets  of  points  should  not  be  expected  to  fall  on  the  same 
Boltzmann  curve.  Obviously,  more  careful  attention  must  be  paid  to  bandpass  (as  in  OH). 
However,  a  small  bandpass  is  required  to  enhance  the  light  rejection  for  propellant 
experiments.  A  better  bandpass  would  be  0.76  nm  (387.58  to  388.34  nm),  which  includes  all 
the  P  branch  lines  from  N=1  to  56  and  none  of  the  R  branch  lines.  The  CN  LIF  intensity  vs 
bandpass  behavior  is  shown  in  Figure  24. 
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Table  2.  Equations  used  for  CN  violet  system  calculations  (B2Z-X2Z). 


T  =  Te  +  G(v)  +  F(N)  =  Te  +  o>e  (v+0.5)  -  C0cXe(v+0.5)2  +  ByNCN+l^DyN^N+l)2 
where  N  =  K  in  Herzberg  [14] 

Fi(N)  =  BVN(N+1)  -Dyl^N+l)2  +  0.5yN  J=N+0.5 

F2(N)  =  BVN(N+1)  -DyN^N+l)2  -  0.5y(N+l)  J=N-0.5 

then,  energy  (wavenumber)  of  transition  is  v  =  T-T",  for  P  branch  N’-N’^l,  R  branch  NP- 
N"=l 

line  strength:  (Honl-London  factor)  Sp=N,  Sr=N+l  for  A=0  (Z-Z  transition) 

Data  from  Engleman  [16]  (Bv,  Dy,  y)  and  Huber  and  Herzberg  [21]  (cue,  o)exc,  Te): 

For  CN  0,1  band  v'=0,  v"=l: 

Bv’=1.95874  BVA=1.87364  Tc=25751.8  ctn-1 

Dy^ti.SSXlO*6  Dv  "=6.42X10*^ 

V=15.65X10-3  y”=5.96X10*3 

(Oe-2163.90  £^"=2068.59 

aicxe’=20.2  Wcxc'^13.087 

For  0,0  band  (v"=0): 

Bv"*=1.8910 

Dv"=6.41X10-<> 

Y"=6.16X10-3 


Figure  6.  CN  LIF  schematic  for  0,1  band  excitation-0,0  band  detection. 
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Table  3:  Calculated  lines  for  CN  0,1  band  including  Jevons  data  [15].  Doublets  not  resolved 

in  [15]. 


N" 

Pi 

CL 

Ri 

r2 

Jevons  P 

Jevons  R 

0 

420.7711 

420.7715 

420.7691 

i 

420.9071 

420.9071 

420.6987 

420.6992 

420.9040 

420.6957 

2 

420.9704 

420.9706 

420.6233 

420.6239 

420.9670 

420.6212 

3 

421.0307 

421.0311 

420.5448 

420.5457 

421.0270 

420.5423 

4 

421.0880 

421.0885 

420.4634 

420.4645 

421.0850 

420.4613 

5 

421.1423 

421.1429 

420.3791 

420.3803 

421.1400 

420.5542 

6 

421.1935 

421.1944 

420.2917 

420.2931 

421.1910 

420.2909 

7 

421.2418 

421.2428 

420.2015 

420.2030 

421.2390 

420.2008 

8 

421.2870 

421.2881 

420.1083 

420.1100 

421.2860 

420.1065 

9 

421.3291 

421.3305 

420.0121 

420.0140 

421.3270 

420.0109 

10 

421.3683 

421.3698 

419.9131 

419.9151 

421.3660 

419.9108 

11 

421.4044 

421.4060 

419.8111 

419.8133 

421.4020 

<19.8096 

12 

421.4374 

421.4393 

419.7062 

419.7086 

421.4360 

419.7044 

13 

421.4674 

421.4695 

419.5983 

419.6009 

421.4660 

419.5917 

14 

421.4943 

421.4966 

419.4876 

419.4904 

421.4920 

419.4876 

15 

421.S182 

421.5206 

419.3741 

419.3770 

421.5170 

419.3730 

16 

421.5390 

421.5416 

419.2576 

419.2607 

421.5400 

419.2570 

17 

421.5568 

421.5596 

419.1382 

419.1414 

421.5570 

419.1366 

18 

421.5715 

421.5744 

419.0160 

419.0194 

421.5720 

419.0157 

19 

421.5832 

421.5862 

418.8910 

418.8945 

421.5840 

418.8902 

20 

421.5917 

421.5950 

418.7630 

418.7668 

21 

421.5972 

421.6007 

418.6323 

418.6362 

22 

421.5997 

421.6033 

418.4987 

418.5028 

23 

421.5990 

421.6028 

418.3624 

418.3666 

24 

421.5953 

421.5992 

418.2232 

418.2276 

25 

421.5885 

421.5927 

418.0811 

418.0858 

26 

421.5786 

421.5829 

417.9364 

417.9412 

27 

421.5657 

421.5702 

417.7888 

417.7938 

28 

421.5497 

421.5544 

417.6385 

417.6436 

29 

421.5306 

421.5354 

417.4855 

417.4908 

30 

421.5085 

421.5134 

417.3297 

417.3351 

31 

421.4832 

421.4884 

417.1711 

417.1768 

32 

421.4549 

421.4602 

417.0099 

417.0157 

33 

421.4235 

421.4290 

416.8459 

416.8518 

34 

421.3890 

421.3947 

416.6793 

416.6853 

35 

421.3515 

421.3573 

416.5099 

416.5162 

36 

421.3109 

421.3169 

416.3379 

416.3443 

37 

421.2672 

421.2733 

416.1632 

416.1698 

38 

421.2205 

421.2268 

415.9859 

415.9926 

39  • 

421.1706 

421.1771 

415.8059 

415.8128 

40 

421.1177 

421.1244 

415.6233 

415.6303 

41 

421.0618 

421.0686 

415.4381 

415.4453 

42 

421.0028 

421.0098 

415.2503 

415.2577 

43 

420.940V 

420.9479 

415.0598 

415.0674 

44 

420.8756 

420.8830 

414.8669 

414.8746 

45 

420.8074 

420.8150 

414.6714 

414.6792 

46 

420.7362 

420.7439 

414.4732 

414.4813 

47 

420.6619 

420.6698 

414.2726 

414.2808 

48 

420.5847 

420.5927 

414.0695 

414.0778 

49 

420.5043 

420.5125 

413.8639 

413.8723 
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Tabic  3:  Calculated  lines  for  CN  0,1  band  including  Jevons  data  [15]  (cont) 


N" 


so 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 
61 
62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 


Pi 

P2 

Ri 

R2 

Jevons  P 

420.4209 

420.4293 

413.6557 

413.6643 

420.3345 

420.3431 

413.4450 

413.4539 

420.2451 

420.2538 

413.2320 

413.2409 

420.1526 

420.1615 

413.0164 

413.0255 

420.0572 

420.0662 

412.7984 

412.8077 

419.9587 

419.9679 

412.5780 

412.5874 

419.8572 

419.8666 

412.3552 

412.3647 

419.7527 

419.7623 

412.1300 

412.1397 

419.6453 

419.6550 

411.9024 

411.9123 

419.5348 

419.5446 

411.6724 

411.6824 

419.4213 

419.4314 

411.4401 

411.4503 

419.3049 

419.3151 

411.2054 

411.2158 

419.1855 

419.1959 

410.9685 

410.9789 

419.0631 

419.0737 

410.7292 

410.7399 

418.9379 

418.9485 

410.4876 

410.4984 

418.8096 

418.8204 

410.2438 

410.2547 

418.6784 

418.6894 

409.9977 

410.0088 

418.5443 

418.5555 

409.7493 

409.7606 

418.4072 

418.4185 

409.4987 

409.5102 

418.2672 

418.2787 

409.2460 

409.2575 

418.1243 

418.1360 

408.9910 

409.0027 

417.9785 

417.9903 

408.7338 

408.74S6 

417.8297 

417.8417 

408.4745 

408.4865 

417.6782 

417.6903 

408.2130 

408.2252 

417.5237 

417.5360 

407.9494 

407.9617 

417.3664 

417.378 8 

407.6837 

407.6960 

417.2061 

417.2188 

407.4158 

407.4284 

417.0431 

417.0559 

407.1459 

407.1586 

416.8772 

416.8901 

406.8739 

406.8867 

416.7084 

416.7216 

406.5999 

406.6129 

416.5369 

416.5502 

406.3238 

406.3369 

416.3625 

416.3759 

406.0457 

406.0590 

416.1854 

416.1989 

405.7656 

405.7791 

416.0054 

416.0191 

405.4836 

405.4971 

415.8226 

415.8365 

405.1995 

405.2132 

415.6371 

415.6511 

404.9135 

404.9274 

415.4488 

415.4630 

404.6256 

404.6396 

415.2577 

415.2721 

404.3358 

404.3499 

415.0639 

415.0784 

404.0440 

404.0583 

414.8674 

414.8821 

403.7504 

403.7648 

414.6682 

414.6830 

403.4550 

403.4695 

414.4662 

414.4812 

403.1576 

403.1723 

414.2616 

414.2767 

402.8585 

402.8733 

414.0542 

414.0695 

402.5575 

402.5724 

413.8442 

413.8596 

402.2547 

402.2698 

413.6315 

413.6471 

401.9502 

401.9654 

413.4162 

413.4319 

401.6439 

401.6592 

413.1983 

413.2141 

401.3358 

401.3513 

412.9777 

412.9937 

401.0261 

401.0417 

412.7545 

412.7707 

400.7146 

400.7303 
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Table  4:  Calculated  line  list  for  CN  0,0  band,  including  Engleman  data  [16] 


N" 

Pi 

P? 

Ri 

r2 

B5I3EI 

0 

387.4635 

387.4638 

387.462 

387.462 

l 

387.5793 

387.5793 

387.4026 

387.4030 

387.577 

387.577 

387.401 

387.401 

2 

387.6341 

387.6342 

387.3397 

387.3403 

387.632 

387.632 

387.338 

387.338 

3 

387.6868 

387.6871 

387.2747 

387.2755 

387.685 

387.685 

387.273 

387.273  , 

4 

387.7375 

387.7379 

387.2078 

387.2087 

387.736 

387.736 

387.206 

387.206 

5 

387.7861 

387.7866 

387.1389 

387.1399 

387.784 

387.784 

387.137 

387.137 

6 

387.8327 

387.8334 

387.0679 

387.0691 

387.831 

387.831 

387.067 

387.067 

7 

387.8773 

387.8781 

386.9951 

386.9964 

387.876 

387.876 

386.993 

386.993  • 

8 

387.9198 

387.9207 

386.9202 

386.9216 

387.919 

387.919 

386.918 

386.918 

9 

387.9602 

387.9614 

386.8433 

386.8449 

387.959 

387.959 

386.841 

386.841 

10 

387.9987 

387.9999 

386.7644 

386.7662 

387.997 

387.997 

386.763 

386.763 

11 

388.0350 

388.0364 

386.6837 

386.6855 

388.034 

388.034 

386.682 

386.682 

12 

388.0693 

388.0709 

386.6009 

386.6029 

388.068 

388.068 

386.600 

386.600 

13 

388.1016 

388.1033 

386.5162 

386.5183 

388.101 

388.101 

386.514 

386.516 

14 

388.1318 

388.1336 

386.4295 

386.4318 

388.131 

388.131 

386.429 

13 

388.1599 

388.1619 

386.3409 

386.3434 

388.158 

388.160 

386.338 

386.341 

16 

388.1859 

388.1881 

386.2504 

386.2530 

388.186 

388.188 

386.248 

.  386.251 

17 

388.2100 

388.2122 

386.1579 

386.1606 

388.208 

388.211 

386.155 

386.158 

18 

388.2318 

388.2343 

386.0635 

386.0664 

388.230 

388.233 

386.061 

386.064 

19 

388.2517 

388.2543 

385.9672 

385.9702 

380.250 

388.253 

385.965 

385.968 

20 

388.2695 

388.2722 

385.8690 

385.8721 

388.268 

388.271 

385.866 

385.870 

21 

388.2852 

388.2880 

385.7689 

385.7722 

388.284 

388.287 

385.767 

385.770 

22 

388.2988 

388.3018 

385.6668 

385.6702 

388.298 

388.300 

385.665 

385.668 

23 

388.3103 

388.3134 

385.5629 

385.5665 

388.309 

388.312 

385.577 

385.583 

24 

388.3198 

388.3230 

385.4571 

385.4608 

388.318 

388.321 

385.473 

385.477 

25 

388.3271 

388.3306 

385.3495 

385.3533 

388.326 

388.330 

385.348 

385.369 

26 

388.3324 

388.3360 

385.2399 

385.2439 

388.330 

388.335 

385.238 

385.260 

27 

388.3357 

388.3394 

385.1285 

385.1326 

388.335 

388.339 

385.143 

385.149 

28 

388.3368 

388.3406 

385.0152 

385.0194 

388.335 

388.339 

385.031 

385.036 

29 

388.3358 

388.3398 

384.9001 

384.9045 

388.335 

388.339 

384 . 902 

30 

388.3328 

388.3369 

384.7831 

384.7877 

31 

388.3276 

388.3319 

384.6644 

384.6690 

32 

388.3204 

388.3248 

384.5437 

384.5485 

33 

388.3111 

388.3157 

384.4213 

384.4262 

34 

388.2997 

388.3044 

384.2971 

384.3021 

35 

388.2862 

388.2911 

384.1710 

384.1762 

36 

388.2707 

388.2756 

384.0432 

384.0485 

37 

388.2530 

388.2581 

383.9135 

383.9190 

38 

388.2333 

388.2385 

383.7821 

383.7878 

39 

388.2115 

388.2169 

383.6489 

383.6547 

40 

388.1875 

388.1931 

383.5140 

383.5199 

41 

388.1616 

388.1672 

383.3773 

383.3833 

42 

388.1335 

388.1393 

383.2389 

383.2450 

43 

388.1033 

388.1093 

383.0987 

383.1050 

44 

388.0710 

388.0772 

382.9568 

382.9633 

4 

45 

388.0367 

388.0430 

382.8132 

382.8198 

46 

388.0004 

388.0068 

382.6679 

382.6746 

47 

387.9619 

387.9684 

382.5209 

382.5277 

48 

387.9214 

387.9280 

382.3722 

382.3791 

49 

387.8787 

387.8856 

382.2218 

382.2289 

50 

387.8340 

387.8410 

382.0697 

382.0769 

12 


Table  4:  Calculated  line  list  for  CN  0,0  band,  including  Engleman  data  [16]  (corn). 


• 

N" 

Pi 

P2 

Ri  _ 

■K9H 

51 

387.7873 

387.7944 

381.9160 

381.9234 

52 

387.7385 

387.7457 

381.7606 

381.7681 

53 

387.6876 

387.6950 

381.6036 

381.6112 

54 

387.6346 

387.6422 

381.4449 

381.4527 

55 

387.5796 

387.5873 

381.2847 

381.2925 

.  56 

387.5226 

387.5304 

381.1227 

381.1308 

57 

387.4635 

387.4714 

380.9593 

380.9674 

58 

387.4023 

387.4104 

380.7942 

380.8024 

59 

387.3391 

387.3473 

380.6275 

380.6359 

'  60 

387.2738 

387.2822 

380.4592 

380.4678 

61 

387.2065 

387.2151 

380.2894 

380.2981 

62 

387.1372 

387.1459 

380.1180 

380.1269 

63 

387.0658 

387.0746 

379.9451 

379.9541 

64 

386.9924 

387.0014 

379.7707 

379.7798 

65 

386.9171 

386.9261 

379.5948 

379.6039 

66 

386.8396 

386.8488 

379.4173 

379.4266 

67 

386.7602 

386.7695 

379.2383 

379.2477 

68 

386.6787 

386.6882 

379.0578 

379.0674 

69 

386.59S2 

386.6048 

378.8759 

378.8856 

70 

386.5097 

386.5195 

378.6924 

378.7023 

71 

386.4223 

386.4322 

378.5076 

378.5175 

72 

386.3328 

386.3428 

378.3213 

378.3313 

73 

386.2413 

386.2515 

378.1335 

378.1437 

74 

386.1479 

386.1582 

377.9443 

377.9546 

75 

386.0525 

386.0629 

377.7537 

377.7642 

76 

385.9550 

385.9656 

377.5617 

377.5723 

77 

385.8557 

385.8664 

377.3683 

377.3790 

78 

385.7543 

385.7652 

377.1736 

377.1844 

79 

385.6510 

385.6620 

376.9774 

376.9883 

80 

385.5458 

385.5569 

376.7799 

376.7910 

81 

385.4386 

385.4499 

376.5811 

376.5923 

82 

385.3295 

385.3409 

376.3809 

376.3922 

83 

385.2184 

385.2299 

376.1794 

376.1908 

84 

385.1054 

385.1170 

375.9766 

375.9881 

85 

384.9904 

385.0023 

375.7725 

375.7841 

86 

384.8736 

384.8856 

375.5670 

375.5789 

87 

384.7549 

384.7669 

375.3604 

375.3723 

88 

384.6342 

384.6464 

375.1524 

375.1645 

89 

384.5117 

384.5240 

374.9432 

374.9554 

90 

384.3872 

384.3997 

374.7328 

374.7451 

91 

384.2609 

384.2735 

374.5211 

374.5335 

92 

384.1327 

384.1454 

374.3083 

374.3208 

93 

384.0026 

384.0155 

374.0942 

374.1068 

94 

383.8707 

383.8837 

373.8789 

373.8917 

95 

383.7369 

383.7500 

373.6624 

373.6753 

96 

383.6012 

383.6145 

373.4448 

373.4578 

97 

383.4637 

383.4771 

373.2260 

373.2391 

98 

383.3245 

383.3380 

373.0061 

373.0193 

99 

383.1833 

383.1970 

372.7850 

372.7984 

Engle.  Pi 


Engle.  P? 


Engle.  Ri 


Engle.  R? 


383  384  385  386  387  388  389 

Wavelength,  nm 


Figure  7.  CN  0,0  LIF  spectrum  from  0,1  excitation.  0.15  mJ/pulse  in  flame, 
S  1=100  pm  X  1  cm,  S4  =200  pm  (0.09  nm  detection  bandpass),  421.6  nm  excitation.  SI,  S4 
refer  to  front  and  rear  slits  of  double  spectrometer. 
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421.1  421.2  421.3  421.4  421.5  421.6  421.7 

Laser  wavelength,  nm 


Figure  8.  CN  1,0  LIF  excitation  spectrum.  0.6  mJ/pulse  in  flame  (est.),  S  1=100  pm 
X  0.2  cm,  S 4=200  pm  (0.09  nm  detection  bandpass). 
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Laser  energy  In  flame  (approx),  mJ/pulse 
Figure  9.  Laser  power  vs  CN  LIF  intensity. 
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Figure  10.  Comparison  of  model  and  Jevons  data  [15]. 


388.40 


388.20 


388.00 


387.80 


387.60 


387.40 


•P  2 

-Engleman  PI 
■  Engleman  P2 

\ 


■R1  j? 
-R2  ^ 
■Engleman  R1 
■Engleman  R2 


387.5 


386.5  E 
c 


386  a 


385.5  > 

CO 

5 


384.5 


0  5  10  15  20  25  30 

N" 

Figure  11.  Comparison  of  model  and  Engleman  data  [16]. 
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Figure  12.  Comparison  of  model  and  experimental  CN  spectra  #1.  Dye  laser 
bandwidth  is  0.2  cm*1  (0.004  nm),  same  as  model.  Model  yields  Boltzmann  distribution,  so 
intensities  aren't  expected  to  match  experimental  LIF  data. 
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17 


Figure  15.  Boltzmann  plot  #1  unsatu rated-high  power  (0.6mJ/pulse  in  flame  (est.)). 


Figure  16.  Boltzmann  plot  #2  saturated-high  power.  A  fit  through  all  points  gives  a 
temperature  of  5230  K. 
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Figure  17.  Boltzmann  plot  #3  unsaturated-low  power  (0.26  mJ/pulse  in  flame). 
S  1=100  fim,  S 4=500  pm,  0.22  nm  bandpass  detection. 


Figure  18.  Boltzmann  plot  #4  saturated,  low  power.  A  fit  through  all  points  gives  a 
temperature  of  4955  K. 
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Figure  19.  CN  rotational  level  population  as  a  function  of  temperature.  The 
equation  used  is  f;=(hcBv(2N+l)/kT)exp(-hcBvN(N+l)/kT)  [9].  Maximum  is  found  by 
solving  2-  { (2N+ 1  )*(hcB v/kT) }  =0  for  given  T. 


419.8  419.9  420  420.1  420.2  420.3  420.4  420.5 


Laser  wavelength,  nm 


Figure  20.  CN  1,0  LIF  excitation  scan.  Laser  power  varied  from  0.04  mJ/pulse  at 
420  nm  to  0.27  mJ/pulse  at  421.9  nm,  S  1=200  pm  X  0.2  cm,  S4=400  |im  (0.18  nm  detection 
bandpass). 
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420.5  420.6  420.7  420.8  420.9  421  421.1  421.2 

Laser  wavelength,  nm 


Figure  21.  Same  conditions  as  Figure  20.  A  third  scan  (in  addition  to  Figures  20  and 
21)  looks  essentially  identical  to  Figure  8  and  is  not  shown. 


419.5  420  420.5  421  421.5 

Model  wavelength,  nm 

Figure  22.  Comparison  of  calculated  and  experimental  CN  0,1  band  line 
locations  (for  experiment  of  Figures  20  and  21). 


21 


C  -10 


-y  -  -7.3292  +  -0.001 1057x  R-  0.9535 
1300  K 


1/90- 3to5 


1000  2000  3000  4000  5000 

Rotational  enargy  [B  N(N+1)J,  cm*1 


Figure  23.  Boltzmann  plot  for  excitation  scans  of  Figures  20  and  21. 
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Figure  24.  CN  LIF  intensity  as  a  function  of  bandpass.  A  possible  explanation  for 
this  behavior  First  linear  region  represents  P  branch  bandhead,  second  linear  region  represents 
slower  inclusion  of  CN  rotational  lines  with  increasing  bandpass.  Similar  behavior  seen  in 
emission,  although  second  linear  region  extends  out  to  larger  detection  bandpasses  because  of 
the  larger  P  bandhead  in  emission  (wider  rotational  distribution). 
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C2  LIF  MEASUREMENTS 


During  the  initial  stages  of  attempting  NO  LIF  at  about  255  nm,  C2  LIF  experiments 
were  attempted  at  about  510  nm.  Coumaxin  334  dye  was  used  with  a  nominal  dye  laser  output 
of  506-537  nm.  It  was  not  too  difficult  to  obtain  C2  LEF  in  the  reaction  zone  of  a  CH4/N2O 
(about  1:2)  flame,  and  even  stronger  C2  LIF  signals  were  seen  in  a  similar  C3H8/N2O  flame. 
Initially,  R  and  P  branch  transitions  from  512  to  514.5  nm  were  excited  and  the  LIF  observed 
in  the  0,0  bandhead  at  516.5  nm.  This  is  illustrated  in  Figures  25-31  for  the  C3H8/N2O  flame 
and  Figure  32  for  the  CH4/N2O  flame.  By  comparing  Figures  29  and  32,  the  lower  C2  LIF 
signal  (concentration)  in  the  CH4  flame  can  be  seen.  The  goal  in  this  work  was  to  obtain 
species  profiles  in  propellant  flames.  An  excite  0,0-observe  0,0  LIF  scheme  will  not  work  in 
propellant  flames  (too  hard  to  discriminate  between  large  scattering  at  laser  wavelength  and 
desired  LIF  signal),  so  one  possibility  is  excite  0,0-observe  0,1  scheme  (the  reverse  of  the 
propellant  CN  LIF  experiments,  which  were  excite  0,1 -observe  0,0).  A  part  of  the  C2  0,1  LIF 
spectrum  from  0,0  excitation  is  shown  in  Figure  33.  An  excitation  scan  is  shown  in  Figures 
34  and  35  for  0,1  detection.  Note  the  significantly  lower  signal  due  to  the  lower  transition 
probability  for  0,1  vs  0,0  LIF.  Detecting  the  LIF  in  the  0,1  bandhead  does  allow  the  excitation 
of  the  0,0  bandhead,  as  shown  in  Figure  35.  For  comparison  purposes,  the  C2  emission  - 
spectrum  in  a  CH4/N2O  flame  is  shown  in  Figure  36. 

C2  has  a  small  rotational  constant  (Bv"=1.7495  cm'1)  so  many  rotational  levels  are 
populated.  The  R  branch  consists  of  triplets  and  the  P  branch  of  doublets,  with  the  splitting 
decreasing  with  increasing  J.  The  rotational  level  with  the  maximum  population  is  shown  in 
Figure  37  as  a  function  of  temperature.  Also  shown  in  Figure  37  is  the  rotational  level  with  a 
population  with  the  minimum  temperature  dependence.  The  equations  used  in  these 
calculations  are  listed  in  Table  5.  At  flame  temperatures  of  2000-2500  K,  the  best  transition  to 
excite  for  C2  profile  measurements  is  about  J=30.  C2  line  locations  are  listed  in  Table  6  (data 
from  Shea  [19]).  I  also  calculated  a  (rough)  rotational  temperature  from  the  excitation  scans  of 
Figures  26-32.  I  assumed  unsaturated  LIF,  and  found  T  from  plotting  In  LIF/Sj*g  vs  the 
rotational  energy  (BVJ(J+1)).  The  temperatures  ranged  from  1300-1600  K,  depending  upon 
the  assumptions  made  for  calculating  the  P  branch  peak  heights  (the  P  branch  peaks  are  small 
and  partially  (at  best)  resolved).  An  example  is  shown  in  Figure  38,  with  the  assumption  that 
the  P  branch  peaks  are  totally  resolved.  This  temperature  is  rather  lower  than  the  expected 
2500  K  (from  OH  LIF  measurements),  but  is  close  to  the  CN  rotational  temperature 
measurements  discussed  in  the  previous  section.  More  research  would  have  to  be  done  to  get 
reliable  C2  temperature  measurements  although  relative  C2  concentration  profiles  can  be 
obtained  with  this  information.  Note  that  the  rotational  temperature  of  the  upper  state  as  shown 
in  the  emission  spectrum  is  much  higher  than  that  of  the  LIF  excitation  scan.  This  is  typical 
behavior  for  chemiluminescence  [20]  where  emission  temperatures  of  5000  K  are  not 
uncommon  for  radicals  such  as  CN  and  C2. 

For  the  propellant  experiments,  the  detection  of  C2  0,1  LIF  with  0,0  bandhead 
excitation  (516.52  nm)  was  attempted  in  a  1  atm  AP  propellant  flame.  The  AP  propellant  was 
chosen  because  the  binder  is  polybutadiene,  and  if  any  propellant  flame  is  going  to  have  C2  in 
detectable  quantities,  it  should  be  the  AP  propellant  However,  when  the  spectrometer  slits 
were  closed  enough  to  discriminate  between  the  scattering  off  the  side  of  the  propellant  and  the 
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desired  C2  LIF  signal,  no  C2  LIF  signal  was  seen.  This  was  not  unexpected  since  C2  was  not 
seen  in  emission  in  propellant  flames  [5],  but  the  effort  was  worthwhile  since  C2  is  probably 
indicative  only  of  binder  combustion.  The  source  of  other  molecules,  such  as  CN,  is  more 
ambiguous. 


Table  5.  General  spectroscopic  information  for  C2  Swan  bands  [21]. 

Transition:  A  3ng  -  X  3Tlu  [22],  d  3Tlg  -  a  ^  [21],  A'=A"=1 

Line  strength  (Honl-London  factors)  (Sj)  [14]: 

S  jR=(J"+ 1  +A")(J"+ 1  -  A")/(  J"+ 1 ) 

SjP*(J"+A")(J"-A")/J" 
g  (degeneracy)=s2J+l 

Rotational  peculation  [18]: 

fj=(hcB  v/kT)(2J+ 1  )(exp-  [B  v  J( J+ 1  )hc/kT] ) 

j*2+j*.(k/hcBv)TaVB=0  (J*=rotational  level  with  minimum  T  dependence) 
Jmax=((2kT/Bvhc)0-M)/2  (Jmax==r’°tational  level  with  maximum  population) 


Table  6.  Shea  data  for  C2  0,0  line  locations  [19]  from  511.5-516.5  nm. 


Ri 

r2 

R3 

Pi 

_  P2 

J" 

515.6117 

1 

2 

3 

515.5024 

4 

515.3783 

5 

515.2972 

515.2535 

6 

515.1531 

515.1202 

7 

515.0175 

514.9808 

514.9352 

8 

514.8636 

514.8350 

514.7952 

9 

514.7150 

514.6837 

514.6487 

10 

514.5507 

514.5264 

514.4951 

11 

514.3892 

514.3622 

514.3362 

12 

514.2135 

514.1923 

514.1678 

13 

514.0411 

514.0173 

513.9978 

14 

513.8538 

513.8340 

513.8139 

516.5292 

516.5116 

15 

513.6682 

513.6468 

513.6294 

516.5292 

516.5116 

16 

513.4698 

513.4510 

513.4342 

516.5116 

516.4961 

17 

513.2726 

513.2523 

513.2386 

516.4961 

516.4814 

18 

513.0661 

513.0445 

513.0295 

516.4814 

516.4556 

19 

512.8516 

512.8337 

512.8189 

516.4425 

516.4268 

20 

24 


Table  6  (corn.).  Shea  data  for  C2  0,0  line  locations  [19]  from  51 1.5-516.5  nm. 


Ri 

R2 

R3 

Pi 

P2 

J” 

512.6294 

512.6141 

512.6005 

516.3870 

21 

512.4079 

512.3911 

512.3821 

516.3606 

516.3464 

22 

512.1732 

512.1577 

512.1469 

516.3096 

516.2939 

23 

511.9416 

511.9237 

511.9193 

516.2563 

516.2404 

24 

511.6935 

511.6790 

511.6699 

516.1878 

516.1735 

25 

511.4500 

511.4325 

511.4325 

516.1219 

516.1084 

26 

516.0417 

516.0282 

27 

515.9691 

515.9502 

28 

515.8695 

515.8562 

29 

515.7783 

515.7652 

30 

515.6715 

515.6591 

31 

515.5678 

515.5561 

32 

515.4476 

515.4357 

33 

515.3323 

515.3199 

34  . 

515.1974 

515.1871 

35 

515.0703 

515.0589 

36 

514.9236 

514.9119 

37 

514.7847 

514.7729 

38 

514.6235 

514.6121 

39 

514.4720 

514.4609 

40 

514.2976 

514.2870 

41 

514.1339 

514.1234 

42 

513.9481 

513.9360 

43 

513.7719 

513.7611 

44 

513.5718 

513.5610 

45 

513.3839 

513.3746 

46 

513.1722 

513.1614 

47 

512.9726 

512.9621 

48 

49 

512.5287 

512.5287 

50 

512.2868 

512.2868 

51 

512.0743 

512.0674 

52 

511.8203 

511.8123 

53 

511.5900 

511.5824 

54 

25 


S15  5 15.4  515.8  516.2  516.6  517 

Wavelength,  nm  (approx) 

Figure  25.  C2  LIF  spectrum  from  0,0  Ri  17  excitation  at  513.47  nm.  Detection 
bandwidth=0.24  nm,  slit=600  |im  X  2  cm,  3  mJ/pulse  @  10  Hz,  10  pulse  avg,  30  ns  gate 
width,  50  mV  sensitivity  on  SR  250.  Conditions  same  for  Fig.  26-35  unless  otherwise  noted. 


512.1  512.2  512.3  512.4  512.5 

Laser  wavelength,  nm 

Figure  26.  C2  LIF  excitation  scan.  Detection  «516.4±0. 12  nm. 
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Laser  wavelength,  nm 


Figure  32.  C2  LIF  excitation  scan,  CH4/N2O  flame. 


Laser  wavelength,  nm 

Figure  35.  C2  LIF  excitation  scan,  detection  »563.4±0.12  nm,  10  mJ/pulse. 
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In  {LIF/S,*(2J+1)}  J  with  maximum  population 


500  1000  1500  2000  2500  3000  3500  4000 

Temperature,  K 


Figure  37.  C2  rotational  level  calculations. 


Figure  38.  C2  rotational  temperature  calculations. 


NO  LEF  MEASUREMENTS 


Current  SHG  crystal  limitations  in  the  FL2002  limits  UV  dye  laser  output  to 
wavelengths  greater  than  270  nm.  The  most  likely  candidate  for  NO  LIF  excitation  within  this 
limitation  is  the  NO  0,4  band  (Table  7,  Figure  39).  NO  LEF  from  the  NO  0,4  band  was 
observed  in  a  CH4/N2O/NO  flame,  although  the  LIF  signal  was  relatively  weak  (-1000  X 
weaker  than  OH  1,1  LIF  observed  via  OH  1,0  excitation).  NO  LIF  was  observed  in  the  0,2 
band  (246±2.25  nm)  while  probing  the  0,4  band  from  271-272  nm.  Detecting  NO  LIF  in  the 
0,0  band  at  226  nm  was  not  very  successful,  probably  due  to  detector  response  limitations.  An 
example  of  the  NO  LEF  spectrum  resulting  from  excitation  of  a  NO  0,4  transition  is  shown  in 
Figure  41.  Examples  of  NO  LEF  excitation  scans  are  shown  in  Figures  42-48.  Excitation  scan 
lines  were  identified  by  comparison  with  synthetic  spectra  created  from  the  equations  in  Tables 
8  and  9.  The  synthetic  spectra  were  created  by  scanning  a  delta  function  (the  laser  line)  over 
the  NO  0,4  band  lines  (N=l-80),  with  the  lines  approximated  as  a  Gaussian  function  with  a 
bandwidth  (FWHM)  of  0.4  cm*1  (0.003  nm).  A  more  accurate  representation  would  be  to  scan 
a  Gaussian  laser  beam  (0.1  cm*1  bandwidth)  over  a  Voigt  profile  (combination  of  a  Gaussian 
and  a  Lorentzian  profile),  but  this  wasn't  necessary  for  line  ID  purposes.  The  experimental 
NO  lines  matched  well  with  the  calculated  lines,  as  shown  in  Figure  49.  The  match  wouldbe  - 
even  better  if  more  care  was  taken  in  correcting  the  FL2002  counter  reading.  Sample  calculated 
spectra  are  shown  in  Figures  50-53.  Note  the -interesting  coincidence  that  the  NO  y  bands  have 
the  same  electronic  transition  as  OH  (A^X2!!).  The  rotational  constants  for  NO  are  much 
smaller  than  those  for  OH,  so  the  NO  rotational  energy  levels  are  much  closer  together  and  thus 
more  levels  are  populated  at  a  given  temperature  than  is  the  case  for  OH.  The  nomenclature  for 
the  NO  lines  is  inconsistent,  so  the  lines  were  labeled  by  the  same  criteria  as  OH  This  differed 
from  Reference  26,  with  R12,  021.  Ql2,  P21.  Pl2»  R21  [26]  =  R21,  Q12.  Q21.  Pl2,  O12,  S21 
(this  paper).  The  line  intensity  in  the  model  spectra  was  calculated  as  the  product  of  the  line 
strength  (or  Honl-London  factor.  Table  9)  and  a  temperature-dependent  Boltzmann  factor  (exp- 
{hcG(v)/kT}). 

Once  the  lines  were  identified,  a  Boltzmann  plot  could  be  created  and  an  NO  rotational 
temperature  could  be  calculated.  The  data  used  in  this  calculation  is  given  in  Tables  8  and  9 
and  uses  identified  lines  in  Figures  45-48  (listed  in  Table  10).  Boltzmann  plots  are  shown  in 
Figure  54  (assuming  zero  saturation)  and  in  Figure  55  (assuming  complete  saturation).  As  can 
be  seen  from  Figures  54  and  55,  the  rotational  temperature  of  NO  was  apparently  well  below 
1000  K,  although  the  data  scatter  is  large.  The  temperatures  would  be  better  fit  by  assuming 
partial  saturation  [1 1],  although  the  temperatures  would  still  be  low. 

To  improve  light  rejection  for  propellant  LIF  measurements,  a  Spex  1404  0.85-m 
double  spectrometer  (w/  an  1800  gr/mm  grating)  was  used  to  reproduce  the  measurements. 
The  maximum  bandpass  of  this  spectrometer  was  measured  to  be  -1  nm  (vs  *4.5  for  the  0.5 
m/1200  gr/mm  grating  and  *20  nm  for  the  0.5  m/300  gr/mm  grating).  Thus  the  NO  LIF 
signals  were  correspondingly  lower  than  those  found  with  the  0.5-m  spectrometer,  as  can  be 
seen  by  comparing  Figure  56  (0.85-m  spectrometer)  to  Figure  44  (0.5  m  spectrometer).  The 
NO  LEF  signals  were  so  small  that  NO  LIF  measurements  in  solid  propellant  flames  were  not 
attempted  with  this  technique. 
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An  OH  LIF  excitation  scan  was  performed  at  the  same  flame  location  as  Figure  18. 
The  results  are  shown  in  Figure  57.  Note  that  the  1-nm  bandpass  contributes  to  considerable 
nonlinearity  in  the  Boltzmann  plots.  Also,  the  OH  signal  was  so  large  compared  to  the  NO 
signal  that  detector  saturation  is  a  possibility.  In  any  case,  the  OH  rotational  temperature  is 
substantially  higher  than  the  NO  rotational  temperature.  This  seems  unreasonable,  so  the  low 
NO  rotational  temperatures  must  be  considered  suspect  for  now.  Zabamick  [1 1]  has  observed 
both  NO  0,4  and  NO  0,0  in  a  CH4/NO2/O2  flame  at  50  torr  and  found  the  profiles  to  be 
similar,  but  not  identical 

A  better  scheme  might  be  excite  0,2  -  observe  0,1  or  0,0.  The  relative  population  of  the 
various  levels  is  shown  in  Figure  58,  demonstrating  that  the  concentration  of  NO  molecules  in 
V"=2  can  be  significant  at  higher  temperatures,  with  the  corresponding  increase  in  LIF  signal. 
The  relative  populations  in  die  various  vibrational  levels  can  be  used  to  calculate  a  vibrational 
temperature  in  emission,  as  shown  in  Figure  59.  The  Boltzmann  plot  is  not  very  linear,  but  the 
v"=2-4  bands  form  a  linear  region  with  a  slope  that  corresponds  to  about  5000  K.  Information 
of  this  type  is  used  to  identify  the  reactions  leading  to  the  electronically  excited  species  (see 
next  section).  Note  that  this  type  of  data  analysis  is  not  appropriate  for  Figure  41,  where  a 
laser  is  used  to  populate  v'=0.  In  this  case,  the  intensity  of  the  emitted  bands  should  be  - 
proportional  to  hvA  [18],  where  v  is  the  frequency  of  the  band  and  A  is  the  Einstein  coefficient 
of  spontaneous  emission.  In  other  words,  band  intensity/VA  -  constant.  The  data  for  the 
bands  in  Figure  41  is  shown  in  Table  11  and  Figure  60.  The  band  intensity/vA  is  not  constant 
(-exponential  with  v"  from  Figure  60),  for  reasons  that  are  not  clear  at  this  time. 

Table  7.  Strong  NO  7  bandheads  (O12,  P2)  [22] 


X,  nm 

Intensity 

v*.  v" 

300.88,  299.76 

4 

0,6 

285.95,  284.98 

7 

0,5 

272.22,  271.32 

8 

0,4 

259.57,  258.75 

9 

0,3 

247.87,  247.11 

10 

0,2 

237.02,  236.33 

10 

0,1 

226.94,  226.28 

8 

0,0 

215.49.  214.91 

7 

_ L.Q _ 
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Table  8.  Energy  level  formulas  for  NO  [26] 


T,=G(v")+Fi(v"J") 

F2,1=B(v")[(J"+0.5)2-1±0.5(X(v")2-4X(v”)+4(J,,+0.5)2)0-5] 
X(v”)=A(v”)/B(vM);  S"=0.5,  1.5 

G(v")= 1904.405(v"+0.5)- 14.187  (v"+0.5)2+0.0240(v"+0.5)3 
A(v")= 1 23 .26-0. 1906(v"+0. 5)-0.0 1 08(v"-K).5)2 
B(v")=1.70427-0.01728(v"+0.5)-0.000037(v"-K).5)2 


T=Te+G(v,)+Fi(v',r) 

Fi(v,,J,)=B(v,)(r-0.5)(J’+0.5) 

F2(vy>B(v')(r+o.5)(r+i.5) 

Tc=43906.37 

G(v,)=2374.307(v’+0.5)-16.106(v,+0.5)2-0.04645(v'+0.5)3 
B(v>1. 99478-0.01 8328(v'+0.5) 


Table  9.  Honl-London  (line  strength)  factors  (H)  [26] 


R2<+)  (2J"+1)2±(2J"+1)U(4J"2+4J,,+1-2X) 

R2l(-)  32(J"+1) 

Rl2(-)  (2J>l)2±(2J"+l)U(4r2+4J’’-7+2*.) 

Rl(+)  32(J"+1) 

P2(+)  (2JI!+l)2±(2J',+l)U(4J"2+4J"-7+2X) 

P2l(-)  32J" 

Pl2(-)  (2J"+l)2±(2J"+l)U(4J"2+4J',+l-2l) 

Pl(+)  32J" 

Q2(+)  (2J" +l)r(4J"2+4J”-l)±U(8J"3+12J"2-2J"+L-2&Yl 

02i(-)  32jH(r+i) 

Ql2(-)  (2J"+ 1 )[  (4J,,2+4J"- 1  )±U(8  J"3+ 1 2J"2-2J"-7 +2X)1 

Ql(+)  32J"(J"+1) 

UKXVMXCv'O-KTT+l)2]-^,  X(v")  as  in  Table  8 
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Table  10.  NO  0,4  band  data  (nonresolved  lines:  R2  9  +  P2 19  @  271 .425  nm,  R2  7  +  Q2 1 1  @ 
271 .588  nm,  P2 12  +  R2  3  @  271 .868  nm) 


D 

N 

Line  Strenath 

9" 

g' 

Calc.  X,  nm 

Q2  14 

14 

4.73620 

28 

428.30560 

271.35999 

Oi227 

27 

1.49840 

54 

1347.73999 

271.37961 

Q2  13 

13 

4.32210 

26 

26 

380.56000 

271.43869 

012  26 

26 

1.49740 

52 

50 

1257.34998 

271.46921 

P218 

18 

3.45040 

36 

34 

652.15002 

271.50470 

R28 

8 

1.02460 

16 

18 

191.16000 

271.51041 

Q2  12 

12 

3.91730 

24 

24 

336.10001 

271.51581 

P2  17 

17 

3.22760 

34 

32 

591.26001 

271.57779 

Oi224 

24 

1.48970 

48 

46 

1086.38000 

271.63321 

P2  16 

16 

3.00860 

32 

30 

533.65997 

271.64581 

q2  10 

10 

3.13750 

20 

20 

257.04999 

271.65500 

R2  6 

6 

0.68760 

12 

14 

138.44000 

271.66821 

012  23 

23 

1.48280  ' 

46 

44 

1005.81000 

271.70761 

P2  15 

15 

2.79350 

30 

28 

479.34000 

271.70880  - 

Q29 

9 

2.76330 

18 

18 

222.46001 

271.71710 

r2s 

5 

0.52960 

10 

12 

117.02000 

271.73959 

P2  14 

14 

2.58270 

28 

26 

428.31000 

271.76682 

Ol2  22 

22 

1.47350 

44 

42 

928.52002 

271 .77689 

Q2  8 

8 

2.40010 

16 

16 

191.16000 

271.77429 

R2  4 

0.37940 

8 

10  . 

98.89000 

271.80600 

P2  13 

2.37640 

26 

24 

380.56000 

271.81970 

Q27 

2.04800 

14 

14 

163.14999 

271.82639 

012  21 

1.46170 

42 

40 

854.51001 

271.84119 

Q26 

1.70710 

12 

12 

138.44000 

271.87360 

012  20 

1 .44720 

40 

38 

783.78003 

271.90030 

P2  11 

1.97830 

22 

20 

294.92999 

271.91061 

Q25 

1.05630 

10 

10 

117.02000 

271.91571 

P210? 

10 

1.78710 

20 

18 

257.04999 

271.94849 

Ol2  2 

2 

0.48610 

4 

2 

72.52430 

272.10049 

0-J2  15 

15 

1.35690 

30 

28 

479.33981 

272.12000 

°12  3 

3 

0.56680 

6 

4 

84.05980 

272.13199 

Oi214 

14 

1.28910 

28 

26 

428.30560 

272.14871 

Ol2  4 

4 

0.65510 

8 

6 

98.89060 

272.15869 

012  13 

13 

1.24860 

26 

24 

380.55820 

272.17239 

°125 

5 

0.74160 

10 

8 

117.01630 

272.18030 

°12  12 

12 

1.20350 

24 

22 

336.09851 

272.19101 

°126 

6 

0.82370 

12 

10 

138.43629 

272.19690 

O12II 

11 

1.15350 

22 

20 

294.92761 

272.20459 

©12  7 

m 

0.90060 

14 

12 

163.14999 

272.20850 

0^2  10 

mSm 

1.09840 

20 

18 

257.04639 

272.21310 

Ol2  8 

0.97210 

16 

14 

191.15680 

272.21509 

_ 9 12,0 _ 

1.03800 

18 

16 

222.45590 

272.21661 
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Figure  39.  NO  emission  spectrum.  CH4/N2O/NO  flame,  1  atm  air,  Reticon  detector, 
1000  data  points  (composite  of  2  700  pixel  spectra),  0.07  nm/pixel  dispersion,  500/200  pm 
slits,  0.5  s  exp.  time. 


Figure  40.  Same  as  Figure  39,  corrected  for  Reticon  wavelength  response 
(Figure  76). 
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Detection  wavelength,  nm  (4.5  nm  bandwidth) 


Figure  41.  NO  LIF  spectrum.  0.1  mJ/pulse  in  flame,  271.19  nm  excitation, 
spectrometer  scanned  at  0.05  nm/s,  SI  =200  pm,  S2=3  mm,  10  ns  gate  width.  S2=rear  slit  of 
0.5  m  spectrometer. 
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Figure  42.  NO  LIF  excitation  scan.  0.1  mJ/pulse,  10  ns  gate  width,  246±2.25  nm 
detection,  SI  =200  pm,  S2=3  mm.  Same  conditions  in  Figures  43-48  except  as  noted. 
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Laser  wavelength,  nm 

Figure  43.  NO  LIF  excitation  scan 
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Figure  42.  NO  LIF  excitation  scan. 
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Figure  45.  NO  LIF  excitation  scan. 
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Figure  46.  NO  LIF  excitation  scan. 
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Figure  47.  NO  LIF  excitation  scan. 


Laser  wavelength,  nm 
Figure  48.  NO  LIF  excitation  scan. 


42 


Figure  SO.  Calculated  NO  spectrum  (compare  with  Figure  46). 
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re  51.  Calculated  NO  spectrum  (compare  with  Figure  47). 
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Figure  52.  Calculated  NO  spectrum  (compare  with  Figure  48). 
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Figure  53.  Calculated  NO  spectrum,  same  as  Figure  52  but  Tf=750  K  (vs  3000  in 
Figure  52). 
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Figure  54.  Boltzmann  plot  for  NO  lines  in  Figures  46-49,  assuming  no  saturation. 
Rotational  temperatures:  P2=520  K,  Oi2=680  K 
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Figure  55.  Same  as  Figure  54,  but  saturation  assumed.  Rotational  temperatures: 
P2=1238  K,  Oi2=916  K. 
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Figure  56.  Same  as  Figure  45,  but  0.35  mJ/pulse,  SI  =400  pm,  and  0.85  m 
spectrometer  (1  nm  bandpass). 
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Figure  57.  OH  LIF  excitation  scan  results  from  same  flame  position  as 
Figures  54  and  55.  Same  conditions,  except  0.15  mJ/pulse,  314.8±0.5  nm  detection). 
Rotational  temperatures:  Ri=1515  K,  Ri’=1785  K  (unsat);  Rj=2505  K,  Ri'=865  K  (sat). 


Figure  58.  Boltzmann  distribution  of  vibrational  population  for  NO. 
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Figure  59.  Vibrational  Boltzmann  plot  of  vibrational  band  intensities  in 
emission  from  Figure  40.  The  line  through  the  last  three  points  gives  a  temperature  of 
5000  K.  FC  =  Franck-Condon  factor. 


Table  11.  LIF  calculations  for  Figure  42.  A  -Einstein  coefficient,  v=frequency  of  bandhead. 


v"=0 

v"=l 

v"=2 

v"=3 

LIF  intensity 

40 

180 

280 

315 

detector  sens 

2.2 

3.0 

3.3 

3.5 

v,  cm*1 

44050 

42190 

40320 

38460 

A  (rcl)  [37] 

1.000 

1.372 

1.121 

0.713 

LIF/v*A  (rel) 

1.00 

2.51 

4.55 

7.95 
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HYDRAZINE  DIFFUSION  FLAME  EXPERIMENTS 


INTRODUCTION 

An  ongoing  project  in  the  Combustion  Research  Laboratory  involves  the  development  of  a 
low  pressure  hydrazine  burner  for  the  study  of  hydrazine  combustion  phenomena.  There  are 
several  aspects  of  hydrazine  flames  of  interest.  The  characteristics  of  the  emission  or 
chemiluminescence  from  low  pressure  hydrazine  flames  (especially  the  UV  emission)  is  of  interest 
for  plume  detection  and  tracking.  An  improved  understanding  of  the  physics  and  chemistry  of 
hydrazine/nitrogen  dioxide  diffusion  flames  is  also  of  interest,  with  the  goal  of  improving  the 
performance  of  storable  liquid  propellant  systems.  For  these  reasons,  a  combined 
theoretical/experimental  project  was  undertaken  to  study  hydrazine  combustion. 

This  section  summarizes  some  of  the  experimental  progress  to  date.  Jeff  Amfahr 
constructed  the  hydrazine  burner  and  its  associated  feed  systems  for  fuels  and  oxidizers  and  tested 
the  burner.  He  also  constructed  the  optical  collection  system  for  performing  emission 
spectroscopy  measurements,  and  he  wrote  the  computerized  data  collection  program  for  the 
detector  controller.  This  section  summarizes  results  the  author  collected  with  the  hydrazine  burner. 
The  results  involve  emission  spectroscopy  of  low  pressure  CH4/O2,  CH4/O2/NO2,  and 
CH4/O2/N2H4  diffusion  flames.  Improvements  in  the  hydrazine  feed  system  SSgt  Tracy 
Christensen  suggested  are  also  discussed. 

The  UV  and  visible  emission  in  plumes  is  due  to  electronic  transitions  between  chemically 
produced  excited  states  of  molecules  and  the  ground  states.  The  species  producing  this  emission 
or  chemiluminescence  are  known,  although  the  reactions  producing  the  excited  state  species  and 
the  rates  of  those  reactions  are  much  less  well  known.  One  reason  for  this  is  that  the  reaction 
pathways  producing  the  emission  are  often  not  the  major  pathways  for  the  production  or 
consumption  of  the  species  in  question.  Thus,  the  chemiluminescent  reactions  may  not  be  that 
important  in  the  overall  combustion  behavior  of  the  system,  and  thus  need  not  be  understood  for 
predicting  the  combustion  performance,  for  example.  The  chemiluminescent  emission  from  flames 
has  been  studied  fairly  extensively  [27,28],  and  results  will  be  cited  as  appropriate.  The  emission 
from  the  plumes  of  hydrazine  fueled  rockets  has  also  been  studied  [33].  In  Reference  33,  the  low 
pressure  plume  from  a  monomethylhydrazine/N02  rocket  was  studied,  with  CO,  NO,  OH,  NH, 
and  CN  emission  identified.  The  resolution  was  insufficient  to  determine  rotational  temperatures, 
however.  The  present  work  has  sufficient  resolution  to  calculate  rotational  temperatures. 

EQUIPMENT  AND  PROCEDURES 

The  hydrazine  burner  is  a  concentric  annular  (Burke-Schumann)  diffusion  flame. 
Hydrazine  is  hypergolic  with  oxidizers  such  as  NO2,  so  only  diffusion  flames  can  be  studied.  The 
burner  is  illustrated  in  Figure  61,  the  fuel/oxidizer  feed  system  in  Figure  62,  and  the  overall  system 
in  Figure  63.  Emission  spectroscopy  is  a  line-of-sight  measurement,  with  several  possible 
locations  for  measurements  shown  in  Figure  61  ("Region  1",  etc.).  For  these  preliminary 
measurements,  the  emission  was  collected  from  Region  1,  2  cm  above  the  burner  surface.  The 
inner  and  outer  tube  diameters  of  the  burner  were  1  and  2.2  cm.  The  flame  was  generally  of  the 
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same  shape  for  fuel/oxidizer  flows  near  stoichiometric  or  oxidizer  rich.  Fuel  rich  flames  were 
unusual,  as  shown  in  Figure  61. 

The  flame  is  created  by  igniting  a  CH4/O2  diffusion  flame  at  atmospheric  pressure  and  then 
evacuating  the  chamber  to  the  desired  pressure,  usually  about  50  mm  Hg  absolute  pressure.  Then, 
the  various  other  fiiels/oxidizers  are  gradually  added  to  or  substituted  for  the  CH4  and  the  O2. 
Some  difficulties  were  encountered  with  the  use  of  NO2  and  N2H4,  which  will  be  discussed  in  the 
Results. 

The  flame  emission  was  collected  with  approximately  1:1  magnification  by  a  collection  lens 
system  and  focussed  onto  the  entrance  slit  to  a  Spex  1870  Triplemate  Spectrometer  equipped  with  a 
EG&G/PAR  1420  Reticon  diode  array  with  700  active  elements  (pixels).  The  emission  could  be 
dispersed  onto  the  diode  array  with  either  a  600,  1200,  or  2400  groove/mm  grating  (turret- 
mounted),  with  resulting  dispersions  of  0.07,  0.035,  and  0.0175  nm/pixel  respectively.  This 
system  has  been  extensively  used  in  emission  spectroscopic  studies  of  high  pressure  solid 
propellant  flames  [2,5].  For  these  measurements,  the  entrance  slit  of  the  spectrometer  (which 
defines  the  spatial  resolution  of  the  mission  measurements  in  the  flame)  was  fixed  at  500  pm  wide 
by  1  cm  high.  The  flame  image  is  rotated  by  the  collection  system,  so  the  collection  volume  in  the 
diffusion  flame  is  approximately  1  cm  wide,  500  pm  wide  and  1  cm  in  depth  (the  depth  in  a  line- 
of-sight  measurement  is  the  depth  of  the  flame). 

CH4/O2  FLAME  RESULTS 

The  CH4/O2  diffusion  flame  was  studied  as  a  baseline.  At  1  atm,  the  flame  is  very  white- 
orange  (sooty),  with  the  only  detected  emission  being  that  from  OH.  As  the  pressure  is  lowered, 
the  flame  gradually  begins  to  assume  the  conical  shape  illustrated  in  Figure  61.  The  orange 
disappears  completely  at  pressures  below  about  200  mm  Hg.  At  40  mm  Hg,  the  emission  from 
Region  1  in  the  flame  consists  only  of  OH,  CH,  and  C2  emission.  In  preliminary  measurements  in 
Region  2,  it  appears  that  Region  2  is  a  more  diffuse  version  of  Region  1,  with  the  same  emission 
characteristics.  The  emission  from  Region  1  is  illustrated  in  Figures  64-69.  In  the  figure  captions, 
the  notations  following  the  molecule  (such  as  the  A^X2!!  following  OH)  refer  to  the  transition 
creating  the  emission,  with  the  chemically  excited  (upper)  state  listed  first  and  the  ground  state 
following.  The  pressure  dependence  of  the  flame  structure  is  due  to  the  differing  pressure 
dependencies  of  diffusion  and  reaction  kinetics.  Diffusion  is  inversely  proportional  to  pressure, 
with  the  diffusion  rate  increasing  with  decreasing  pressure.  The  reaction  rate  is  proportional  to  the 
concentration  of  reactants,  and  thus  the  reaction  rate  is  proportional  to  pressure.  Thus,  at 
atmospheric  pressure  where  reaction  is  fast  and  diffusion  is  slow,  the  fuel  is  significantly  reacted 
(pyrolyzed)  before  the  fuel  and  oxidizer  can  mix  through  diffusion.  Hence,  large  amounts  of  soot 
are  formed.  At  low  pressure  diffusion  is  fast,  so  the-fuel  and  oxidizer  can  mix  before  reacting, 
thus  the  reaction  zone  looks  like  a  premixed  flame  with  emission  from  species  such  as  CH  and  C2 
seen  in  the  reaction  zone  of  premixed  flames  [20,27,28].  The  emission  from  Region  1  is 
proportional  to  the  O2/CH4  flow  ratio,  as  shown  in  Figure  70.  This  trend  is  possibly  partially  due 
to  changes  in  flame  shape  and  partially  due  to  other  factors. 
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CH4/O2/NO2  FLAME  RESULTS 


In  premixed  flames,  NO2  is  a  weaker  oxidizer  than  O2  [1 1].  In  the  low  pressure  diffusion 
flame,  it  was  discovered  that  an  CH4/NO2  diffusion  flame  would  not  bum.  When  NO2  was  added 
to  CH4/O2  flame  in  large  excess,  relatively  weak  emission  from  CN  and  NH  indicated  that  a 
relatively  small  amount  of  NO2  was  reacting  in  the  flame.  This  emission  is  illustrated  in  Figures 
71*73.  No  emission  was  seen  in  the  220-280  nm  region,  where  NO  y  band  emission  is  seen  [5]. 
For  NO2  to  react  in  a  low  pressure  diffusion  flame,  a  hypergolic  reaction  such  as  N2H4+NO2  is 
probably  required.  Problems  were  encountered  with  the  NO2  flow  controller,  but  discussions  with 
other  lab  personnel  indicated  that  these  problems  were  due  to  the  low  vapor  pressure  of  NO2  and 
could  be  overcome. 

CH4/O2/N2H4  FLAME  RESULTS 

When  an  attempt  was  made  to  add  N2H4  to  the  CH4/O2  flame,  it  was  found  that  the  feed 
system  illustrated  in  Figure  62  was  inadequate.  The  system  was  designed  to  supply  gaseous  N2H4 
to  the  burner,  but  at  the  conditions  of  the  experiment,  the  N2H4  flowing  through  the  control  valve 
was  liquid  (-298  K,  45  mm  Hg  pressure).  With  liquid  flowing  through  the  lines,  there  was 
negligible  flow  control  of  the  hydrazine.  What  appeared  to  be  happening  was  that  N2H4  was 
entering  the  chamber  as  a  liquid  and  forming  a  pool,  where  it  was  then  flashing  into  a  gas.  As  it 
turned  out,  this  allowed  the  formation  of  a  relatively  stable  CH4/O2/N2H4  flame  of  unknown 
stoichiometry,  illustrated  in  Figure  74.  The  problems  with  the  feed  system  are  being  studied,  with 
one  possible  design  improvement  illustrated  in  Figure  75.  Additional  problems  with  a  leaking 
N2H4  valve  were  also  encountered  and  led  to  an  early  termination  of  the  experiments  with  N2H4. 

The  relatively  stable  CH4/O2/N2H4  flame  was  studied  with  emission  spectroscopy  over  the 
entire  sensitivity  range  of  the  Reticon  (shown  in  Figure  76).  In  addition  to  the  OH,  CH,  and  C2 
emission  from  the  CH4/O2  flame,  strong  CN  and  NH  emission  was  also  seen.  The  emission  from 
the  CH4/O2/N2H4  flame  is  illustrated  in  Figures  77-82.  No  emission  was  detected  from  220-280 
run,  although  several  attempts  with  varying  exposure  times  and  dispersions  were  made.  The  signal 
in  Figure  78  at  220  nm  is  noise  amplified  by  the  detector  sensitivity  correction.  The  apparent 
signal  at  580-600  nm  may  be  NH2  emission  [22],  although  it  may  again  be  amplified  detector 
noise.  The  OH,  CH,  and  C2  emission  are  essentially  unchanged  from  the  CH4/O2  flame,  so  only 
the  CN  and  NH  emission  are  due  to  the  presence  of  the  N2H4.  The  presence  of  NH  in  the  ground 
state  is  reasonable,  with  the  combustion  of  N2H4  proceeding  among  several  competing  paths 
[29,30]: 

N2H4->NH2-*NH->NO->N2  or  N2H4-*N2Hr>NNH->N2 
The  presence  of  CN  indicates  that  the  two  fuels  are  interacting,  since  the  C  comes  from  the  CH4 
and  the  N  from  the  N2H4.  From  the  structure  of  the  flame  as  shown  in  Figure  74,  it  appears  that 
the  reactions  involving  CN  are  occurring  after  the  CH/C2  reactions.  One  possibility  is  that  CH  and 
C2  from  the  CH4/O2  reaction  zone  are  forming  excited  state  CN  and  NH  through  reactions  such  as 
those  with  NO  and  OH  [28]: 


C2(a3II)  +  NO-*CN(B2E)  +  CO 

R1 

CH  +  NO->CN(A2n)  +  OH 

R2 

CN(A2II)  +  OH-*NH(A3ri)  +  CO 

R3 

CH4/N2O  FLAME  RESULTS 

For  comparison,  an  emission  survey  was  made  of  a  CH4/N2O  diffusion  flame.  There  is 
existing  emission  data  taken  with  the  Reticon  detection  system  for  a  premixed,  1  atm  CH4/N2O 
flame  [5],  showing  emission  from  OH,  NH,  CN,  CH,  and  C2  (weak).  The  premixed  flame 
emission  is  illustrated  in  Figure  83  [5].  No  NO  emission  was  seen  in  this  flame.  NO  emission 
was  seen  in  an  NO- seeded  CH4/N2O  flame,  as  shown  in  Figures  40  and  41. 

The  diffusion  flame  of  CH4/N2O  has  slighdy  different  burning  characteristics  from  the 
CH4/O2  diffusion  flame.  The  CH4/N2O  diffusion  flame  was  unstable  below  about  50-60  mm  Hg 
pressure,  so  the  emission  data  was  collected  at  80  mm  Hg  pressure.  The  stoichiometry  of  the 
CH4/N2O  flame  is  different  from  that  of  the  CH4/O2  flame,  with  a  diffusion  flame  similar  to  Figure 
61  (left  drawing)  occurring  at  flow  rates  of  0.4/2.5  slpm  for  CH4/N2O.  The  emission  from 
Region  1  of  this  flame  is  illustrated  in  Figures  84  and  85.  Note  the  similarity  between  the  CN 
emission  shown  in  Figures  85  and  81.  Comparing  the  CH4/N2O  premixed  and  diffusion  flames, 
OH  emission  is  relatively  stronger  in  the  pre mixed  flame,  while  CH  and  C2  emission  is  relatively 
more  intense  in  the  diffusion  flame. 

ROTATIONAL  AND  VIBRATIONAL  TEMPERATURES 

The  rotational  and  vibrational  temperatures  of  the  CN,  NH,  OH,  CH,  and  C2  can  be 
calculated  using  relatively  standard  procedures  previously  employed  in  the  Combustion  Research 
Laboratory  for  analysis  of  solid  propellant  emission  spectra  [5].  These  temperatures  may  give 
some  clues  as  to  the  reactions  producing  the  electronically  excited  CN  and  NR  For  this  report,  the 
main  species  of  interest  are  the  CN  and  NH  present  because  of  the  hydrazine.  The 
rotational/vibrational  temperature  of  experimental  spectra  can  be  calculated  by  comparison  with 
synthetic  spectra,  an  example  of  which  is  shown  in  Figure  86b  [5].  In  Figure  87,  a  corrected  CN 
emission  spectrum  from  the  CH4/O2/N2H4  diffusion  flame  is  shown.  Examples  of  atmospheric 
pressure  premixed  CH4/N2O  flame  CN  spectra  and  high  pressure  propellant  flame  CN  spectra  are 
shown  in  Figures  86a  and  88.  By  comparison  with  synthetic  spectra,  the  CN  vibrational  and 
rotational  temperatures  (assumed  to  be  equal)  for  1  atm  flames  (both  CH4/N2O  and  AP  propellant) 
were  found  to  be  about  5000  K.  The  high  pressure  HMX  propellant  flame  CN  spectra  yielded 
vibrational/rotational  temperatures  of  about  2500  K  [5]. 

By  comparison  with  synthetic  spectra,  the  vibrational/rotational  temperatures  in  the 
CH4/O2/N2H4  diffusion  flame  were  approximately  10500/7000  K  at  45  mm  Hg  pressure  (0.06 
atm).  The  synthetic  spectrum  is  shown  in  Figure  89.  It  was  necessary  to  choose  different 
rotational  and  vibrational  temperatures  to  fit  the  experimental  data  (Figure  87).  These  values  of 
rotational  and  vibrational  temperatures  for  CN(B2£)  are  very  close  to  those  fot  in  a  study  of  the 
C2(a3IT)  +  NO->CN(B2I)  +  CO  reaction  by  Reisler,  Mangir,  and  Wittig  [31]  and  the  rotational 
temperature  of  CN(B2I)  found  in  C2H2/N2O  flames  by  Guillaume  and  Van  Tiggelen  [28] .  Thus, 
the  diffusion  flame  data  at  low  pressures  appears  to  indicate  that  the  CN  emission  in  these  flames  is 
due  to  the  C2+NO  reaction,  as  it  was  in  C2H2/N2 O  flames  [28]. 

The  NH  vibrational  temperature  in  the  CH4/O2/N2H4  diffusion  flame  at  45  mm  Hg  (Figure 
79)  is  about  2750  K,  similar  to  the  temperature  found  in  CH4/N2O  premixed  flames  at  1  atm  and 
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propellant  flames  at  pressures  from  1-35  atm  [5].  A  sample  NH  synthetic  spectrum  is  shown  in 
Figure  90,  with  the  vibrational  temperature  calculated  from  the  ratio  of  the  0,0  and  1,1  Q  branch 
bandheads. 

With  the  addition  of  this  low  pressure  diffusion  flame  data,  a  database  of  CN  and  NH 
vibrational  temperatures  now  exists  over  a  wide  range  of  pressures.  When  the  CN  vibrational 
temperature  is  plotted  as  a  function  of  pressures  for  these  different  flames,  a  strong  correlation  is 
seen,  as  shown  in  Figure  91.  This  correlation  probably  indicates  that  the  reactions  forming  the 
excited-state  CN  (CN*)  are  similar  in  all  these  flames,  with  the  temperature  difference  due  to 
pressure.  What  is  apparently  happening  is  that  the  reactions  are  forming  CN*  in  a  highly- 
vibrationally-excited  state.  At  low  pressures  with  relatively  few  collisions,  this  excited  state 
radiates  with  a  temperature  similar  to  that  at  which  it  was  produced.  At  high  pressures,  the  CN*  is 
involved  in  many  collisions  before  radiating,  with  these  collisions  tending  to  "thermalize"  or  reduce 
the  temperature  of  the  CN*  to  the  flame  temperature  (die  adiabatic  flame  temperature  of  all  these 
flames  is  in  the  2500-3000  K  range).  This  seems  to  also  explain  the  variation  of  CN*  temperature 
with  height  seen  in  AP  propellant  flames  [5],  where  the  temperature  fell  as  distance  above  the 
surface  (and  thus  time)  increased.  This  trend  is  not  seen  with  NH,  mainly  because  NH  is  not 
produced  with  high  vibrational  excitation.  This  is  shown  in  Figure  92. 

The  high  vibrational/rotational  excitation  of  the  CN  and  the  more  "thermal"  distribution  of 
the  NH  is  due  to  the  energetics  of  the  reactions  forming  the  two  species.  The  amount  of  energy 
liberated  or  consumed  in  a  reaction  determines  the  final  state  of  the  species  created.  Using  the  heat 
of  formation  and  excitation  data  from  Table  9,  the  energetics  of  reactions  R1  and  R3  have  been 
calculated: 

C2(a3Il)  +  NO-*CN(B2Z)  +  CO  AH=-74.5  kcal/mole  (exothermic) 

CN(A2n)  +  OH-*NH(A3II)  +  CO  AH=12.2  kcal/mole 

The  heat  of  reaction  of  R1  has  been  calculated  as  -70.5  kcal/mole  by  Guillaume  and  Van  Tiggelen 
[28],  fairly  good  agreement  considering  the  uncertainties  in  the  heats  of  formation  of  the  radical 
species.  Note  that  R1  produces  excess  energy  over  and  above  the  energy  required  to  create  the 
excited  state  CN.  This  excess  energy  appears  in  the  products  as  excess  rotational  and  vibrational 
energy.  The  energy  difference  between  vibrational  levels  in  CN(B2E)  is  approximately  2000  cm'1, 
so  70-75  kcal/mole  (equivalent  to  24,500-26,250  cm*1)  is  easily  enough  energy  to  explain  the 
vibrational  excitation  observed  in  Figures  81  and  86.  In  contrast,  R3  does  not  liberate  excess 
energy  in  forming  NH(A3II),  so  that  the  "thermal"  vibrational  distribution  of  NH  is  consistent  with 
R3  being  the  main  channel  forming  NH(A3II). 
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BURKE-SCHUMANN  CALCULATIONS 


The  "classic"  analysis  of  diffusion  flame  shape  is  that  of  Burke  and  Schumann  [35].  In 
their  analysis  of  diffusion  flames  in  concentric  tubes,  they  obtained  the  following  result  for  over¬ 
ventilated  (fuel-lean)  flames: 


2  {((fci)-1  {Ji(4>iL)  Jo(<M)/[Jo(<l>nR)]2}  exp  [-ItyJyM  1=  (R2C2/2LiCo)-0.5L  (1) 


where  Jo,  Ji=  Bessel  functions 
R=radius  of  outer  tube 
L-radius  of  inner  tube 
x=radial  distance  of  flame  sheet 
v=gas  velocity 

y=veitical  distance  above  orifice  of  inner  tube 
D=fuel-oxidizer  diffusion  coefficient 
Ci=inidal  concentration  of  fuel 
C2=inital  concentration  of  oxidizer 
Co=Ci+(C2fi) 

i=Stoichiometry  coefficient,  moles  oxidizer/moles  fuel 
4>n=root  of  Ji($)=0  (<j>i=3.83) 

The  flame  height  is  calculated  by  setting  x=0  and  solving  for  y  in  Equation  (1).  A  rough  estimate 
can  be  obtained  using  only  n=l  in  the  summation  [35]. 

Data  for  CH4/O2  flame: 

L=0.5  cm 
R=l.l  cm 

Ci=C2=l  (no  diluent) 
i=2  (CH4+2O2— >C02+2H20) 

CH4=0.8  slpm=251.7  cm3/s  @  40  torr=80.1  cm/s 

O2  flow=2.5  slpm=786.3  cm3/s  @  40  torr=65.2  cm/s 

D=4.2  cm2/s  @  40  toir,  298  K  [36],  with  D/v~pD«constant  in  flame 

The  end  result  is  y=7.9  cm,  compared  with  experimental  value  of  3-3.5  cm.  The  agreement  is  not 
as  good  as  that  reported  by  Burke  and  Schumann  for  1  atm  flames  [35].  However,  some 
indication  has  recently  surfaced  that  the  flow  meters  were  reading  inaccurately,  so  this  calculation 
needs  to  be  repeated  with  carefully  calibrated  flow  meters. 
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Table  12.  Energetics  data  for  species  involved  in  reactions  R1 ,  R2,  and  R3. 


Heat  of  formation. 

Heat  of  formation. 

Molecule 

kcal/mole  T321 

kcal/mole  T341 

CO 

-26.42 

-29.4 

NO 

21.6 

20.8 

c2 

199 

198.4 

CN 

101 

100.5 

OH 

9.4 

9.5 

NH 

90 

86.9 

CH 

142 

140.9 

Excited  electronic  state  energies  [21]*: 


Molecule 

State 

Energy,  cm' 

c2 

d^n. 

20022.5 

c2 

aSn,, 

7^.2 

o> 

X’V 

0 

NH 

A3H 

29807 

NH 

X3s- 

0 

CO 

a3nr 

48686.7 

CO 

X  lZ+ 

0 

CN 

B2I+ 

25752 

CN 

A2rii 

-  9245.3 

CN 

X2I+ 

0 

OH 

a2i+ 

32684.1 

OH 

x2nj 

0 

NO 

a2i+ 

43965.7 

NO 

X2rir 

0 

CH 

B2z- 

(26044) 

CH 

A2A 

23189.8 

CH 

X2nr 

0 

n2 

A  3Z+ 

50203.6 

n2 

X  1Z+ 

0 

Common  names  for  transitions:  C2  d-»a  "Swan  bands",  CO  a->X  "Cameron  bands",  CN  A->X 
"CN  red  system",  CN  B-+X  "CN  violet  system",  NO  A-*X  "y  bands" 
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Figure  61.  CH4/O2  diffusion  flame  structure  at  40  mm  Hg  pressure. 
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Figure  62.  Fuel/oxidizer  feed  schematic. 
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Reticon 


Figure  63.  Schematic  (top  view)  of  apparatus  showing  optical  collection  path. 
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Figure  64.  OH  (A2Z-X2FI)  emission  spectrum,  CH4/O2  diffusion  flame  (0.8/2.5  slpm), 
40  mm  Hg  pressure,  0.5  s  exposure  time,  0.0175  nm/pixel  dispersion. 
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Figure  65.  CH  (B2£-X2II)  emission  spectrum,  CH4/O2  diffusion  flame  (0.8/2.5  slpm), 
40  mm  Hg  pressure,  0.2  s  exposure  time,  0.0175  nm/pixel  dispersion. 


Figure  66.  CH  (A2A-X2I1)  emission  spectrum,  CH4/O2  diffusion  flame  (0.8/2.5  slpm), 
0.2  s  exposure  time,  40  mm  Hg  pressure,  0.0175  nm/pixel  dispersion. 
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Figure  67.  C2  (d3Il-a3n)  emission  spectrum,  CH4/O2  diffusion  flame  (0.8/2.5  slpm), 
0.2  s  exposure  time,  40  mm  Hg  pressure,  0.0175  nm/pixel  dispersion. 


Figure  68.  C2  (d3Il-a3n)  emission  spectrum,  CH4/O2  diffusion  flame  (0.8/2.5  slpm),  40 
mm  Hg  pressure,  0.2  s  exposure  time,  0.0175  nm/pixel  dispersion. 
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Figure  69.  C2  (d3n-a3Il)  emission  spectrum,  CH4/O2  diffusion  flame  (0.8/2.5  slpm),  40 
mm  Hg  pressure,  0.2  s  exposure  time,  0.035  nm/pixel  dispersion. 
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Figure  70.  CH  emission  intensity  (A-X  bandhead  at  431  nm)  as  a  function  of  O2 
flow  rate  (CH4  constant  at  0.8  slpm). 
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Figure  71.  NH  (A3I1-X3Z)  emission  spectrum,  CH4/O2/NO2  diffusion  flame  (0.8/2.5/25 
slpm),  40  mm  Hg  pressure,  0.2  s  exposure  time,  0.035  nm/pixel  dispersion. 


Figure  72.  CN  (B2I-X2Z)  emission  spectrum,  CH4/O2/NO2  diffusion  flame  (0.8/2.5/25 
slpm),  40  mm  Hg  pressure,  0.2  s  exposure  time,  0.035  nm/pixel  dispersion. 
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Figure  73.  CN  (B2E-X2X)  (+CH)  emission  spectrum*  CH4/O2/NO2  diffusion  flame 
(0.8/2.5/25  slpm),  40  mm  Hg  pressure,  0.2  s  exposure  time,  0.035  nm/pixel  dispersion. 


Figure  74.  CH4/O2/N2H4  flame  schematic  (45  mm  Hg). 
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Figure  75.  Schematic  of  new  N2H4  feed  system. 
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Figure  76.  Reticon  sensitivity  as  a  function  of  wavelength. 
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Reticon  spectral  response  (Figure  76). 


Figure  79.  NH  (A3n*X3Z)  emission  spectrum,  CH4/N2H4/O2  diffusion  flame  (0.8/7/2.5 
slpm),  45  mm  Hg  pressure,  0.4  s  exposure  time,  0.0175  nm/pixel  dispersion. 


Figure  80.  CN  (B2I-X2Z)  emission  spectrum,  CH4/N2H4/O2  diffusion  flame  (0.8/7/2.5 
slpm),  45  mm  Hg  pressure,  0.4  s  exposure  time,  0.0175  nm/pixel  dispersion. 
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Figure  81.  CN  (B2Z-X2Z)  (+CH)  emission  spectrum,  CH4/N2H4/O2  diffusion  flame 
(0.8/7/2.5  slpm),  45  mm  Hg  pressure,  0.4  s  exposure  time,  0.0175  nm/pixel  dispersion. 


Figure  82.  Emission  spectrum,  CH4/N2H4/O2  diffusion  flame  (0.8/7/2.5  slpm),  50  mm  Hg 
pressure,  0.2  s  exposure  time,  0.035  nm/pixel  dispersion. 
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Wavelength,  nm 

Figure  83.  Broadband  emission  spectrum  of  inner  cone  of  CH4/N2O  premixed 
flame,  1  atm.  500/200  Jim  slits,  0.07  nm/pixel  dispersion,  0.2  s  exposure  time. 


Figure  84.  CN  (B2I-X2I)  emission  spectrum,  CH4/N2O  diffusion  flame  (0.4/2.5  slpm), 
80  mm  Hg  pressure,  0.4  s  exposure  time,  0.0175  nm/pixel  dispersion. 
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Figure  85.  Composite  emission  spectrum.  CH4/N2O  diffusion  flame  (0.4/2.5  slpm),  80 
mm  Hg  pressure,  0.4  s  exposure  time,  0.035  nm/pixel  dispersion. 


Figure  86.  Sample  CN  (B2I-X2Z)  emission  spectra,  (a)  CH4/N2O  flame,  1  atm.  (b) 
Model  spectrum  with  T^T^SOOO  K. 
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Figure  87.  CN  (B2I-X2I)  emission  spectrum,  CH4/N2H4/O2  diffusion  flame,  created  by 
subtracting  CH  emission  (Figure  5)  from  Figure  20. 


Figure  88.  CN  (B2£-X2I)  emission.  HMX1  propellant,  3.5  MPa,  0.0175  nm/pixel 
dispersion,  0.24  s  exp.  time,  500/100  pra  slits,  h=i-2mm. 
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Figure  89.  Model  CN  (B2X-X2X)  emission  spectrum.  Tr=7000  K,  Tv= 10500  K  (0.0 
band  sequence). 


Figure  90.  Model  NH  (A3n-X3X)  emission  spectrum  ^=7^=2500  K). 


Figure  91. 
temperature. 


Figure  92. 
temperature. 


Relationship  between  pressure  and  CN  emission  (upper  state) 


y  -  3401 .4  +  248.42k>g(x)  R-  0.45231 


Relationship  between  pressure  and  NH  emission  (upper  state) 
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CHEMICAL  KINETIC  MODELING 


INTRODUCTION 

Current  solid  propellant  combustion  models  treat  the  propellant  flame  chemistry  in  a 
relatively  limited  way,  defining  the  chemistry  in  terms  of  several  global  reactions,  and 
determining  the  combustion  behavior  mainly  via  a  heat  transfer  analysis  [38,39].  In  some 
cases,  the  lack  of  chemical  detail  has  limited  the  predictive  capability  of  these  type  of  models. 
Conversely,  the  general  lack  of  understanding  of  the  propellant  flame  chemistry  has  limited  the 
development  of  global  reaction  combustion  models,  since  it  is  not  clear  which  global  reactions 
should  be  used.  To  improve  the  knowledge  of  propellant  flame  chemistry  more  research  needs 
to  be  done. 

Validation  of  detailed  models  of  propellant  combustion  requires  comparison  of  model 
predictions  with  experimental  data  collected  in  the  propellant  flame  under  combustion 
conditions  of  interest  The  needed  data  are  the  temperature  profile  and  species  profiles  in  the 
propellant  flame  under  combustion  conditions.  Such  data  is  relatively  sparse  and  difficult  to 
obtain,  compared  with  (for  example)  the  amount  of  information  available  for  methane/air 
flames.  Kubota's  group  in  Japan  has  measured  temperature  profiles  in  propellants  using 
microthermocouples  [40,41].  Other  microthermocouple  measurements  in  propellant  flames 
have  been  performed  in  Italy  [42],  France  [43],  the  Soviet  Union  [44]  and  the  United  States 
[45,64].  Probe  sampling  of  species  has  been  done  in  propellant  flames  in  Japan  [40]  and  the 
Soviet  Union  [46].  Both  microthermocouple  temperature  measurements  and  probe  sampling 
are  intrusive  measurements,  with  the  possibility  of  disturbing  the  flame  chemistry  and  structure 
during  measurement  Nonintrusive  optical  measurements  have  been  made  in  propellant  flames, 
by  such  methods  as  emission  [2],  absorption  [12],  CARS  [13],  and  laser-induced  fluorescence 
(LIF)  [4,5,10].  IR  spectroscopy  has  been  used  to  detect  product  species  in  high  rate 
thermolysis  of  propellant  ingredients  [47].  It  is  this  detailed  data  to  which  model  predictions 
need  to  be  compared,  rather  than  just  comparing  experimental  and  theoretical  bum  rates,  for 
example. 

Detailed  kinetic  modeling  of  propellant  flames  is  not  a  new  development  [48,49], 
although  chemical  kinetic  modeling  in  general  has  only  recently  become  "routine"  with 
production  codes  such  as  the  Sandia  premixed  flame  code  [50]  often  used  to  model  various 
types  of  flames  [30].  Propellant-related  experimental/modeling  efforts  are  underway,  such  as 
CH2O/NO2  and  CH2O/N2O  [51,52],  C2N2/NO2  and  HCN/NO2  [53],  and  CH4/NO2  [11,52]. 
Research  in  the  Soviet  Union  has  also  been  active,  with  kinetic  models  developed  for 
NH3/HCIO4  (AP)  flames  [54]  and  hexogen  (RDX)  flames  [55].  In  the  1960's,  there  was  a 
great  deal  of  research,  primarily  in  Great  Britain  and  France,  on  C10x*oxidized  flames  [56]. 
The  recent  increases  in  computer  capabilities  have  allowed  detailed  kinetic  models  of  nitramine 
combustion  to  be  developed  which  include  the  solid  — »  gas  transition,  thus  allowing  the 
determination  of  the  propellant  bum  rate.  Hatch  developed  complete  models  for  nitrate  ester 
and  HMX  combustion  [57],  the  latter  utilizing  77  reactions  involving  26  species.  Melius 
developed  a  detailed  chemical  kinetic  model  of  RDX  ignition  and  RDX/HMX  combustion  [58], 
with  a  recommended  set  of  131  reactions  involving  32  species  (including  the  nitramine 
molecule).  One  of  Melius'  most  interesting  conclusions  was  that  the  higher  vapor  pressure  of 
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RDX  led  to  very  little  condensed  phase  reaction,  while  over  50%  of  the  higher  boiling  HMX 
reacted  in  the  condensed  phase  prior  to  vaporization.  The  condensed  phase  decomposition  of 
HMX  is  generally  thought  to  produce  CH2O  +  N2O,  while  the  gas  phase  decomposition 
produces  NO2  and  H2CN  through  a  concerted  ring  degradation  mechanism  after  an  NO2  group 
is  removed  from  the  parent  nitramine  [47,59].  This  picture  is  not  universally  agreed  upon, 
with  some  HMX  combustion  models  neglecting  the  H2CN/NO2  path  altogether,  focussing  on 
the  CH2O/N2O  pathway  [60]. 

In  this  section,  available  experimental  species  data  from  solid  propellant  flames  is 
compared  with  chemical  kinetic  models  for  propellants  HMX1,  HMX2,  and  AN1  (see  Table 
1).  HMX2  is  similar  to  a  propellant  studied  by  Kubota  [40],  while  HMX1  has  been  studied  by 
a  variety  of  diagnostic  techniques  by  several  groups  [4,5,10,12,13,61].  AN  propellants  are  of 
current  interest  because  of  the  desire  to  minimize  toxic  exhaust  products,  such  as  the  HC1 
found  in  the  exhaust  of  rockets  burning  ammonium  perchlorate-based  propellants.  Two 
chemical  kinetic  mechanisms  are  used.  The  first  is  the  mechanism  of  Hatch  [57]  with  the 
reaction  H2CN+M~>H+HCN  added.  This  mechanism  is  listed  in  Appendix  1.  The  second 
mechanism  is  the  Miller-Bowman  mechanism  for  nitrogen  chemistry  in  flames  [30],  as 
extended  by  Zabamick  [1 1].  This  mechanism  is  listed  in  Appendix  2.  The  Hatch  mechanism  ' 
doesn't  include  hydrocarbons,  so  binder  decomposition  products  were  limited  to  CH2O.  This 
will  be  discussed  further  in  the  next  subsection.  The  propellants  are  modeled  as  premixed, 
burner-stabilized  flames,  an  approximation  discussed  in  the  following  subsections  This  type  of 
model  requires  a  temperature  profile  as  an  input.  The  codes  were  run  on  a  Macintosh  II 
computer  with  4  megabytes  of  RAM  (Hatch  mechanism)  and  a  VAX  8650  (extended  Miller- 
Bowman  mechanism). 

HMX1  PROPELLANT 

For  the  modeling  portion  of  this  work,  the  decision  was  made  to  narrow  the  focus  of 
the  model  to  the  changes  in  flame  structure  created  by  adding  binder  decomposition  species  to  a 
"baseline"  HMX  reaction  scheme.  To  avoid  solving  the  coupled  gas/solid  problem  with  the 
wide  swings  in  bum  rate  possible  depending  on  the  species  added,  a  simpler  situation  was 
solved.  An  experimental  temperature  profile  was  used  with  the  Sandia  premixed  flame  code  to 
solve  a  burner-stabilized  flame  of  the  same  mass  burning  rate  as  that  of  the  propellant  (which 
was  known  [5]).  Thus  the  temperature  profile  and  bum  rate  are  inputs  (and  assumed  to  be 
correct),  and  the  starting  species  are  varied  to  determine  the  effect  of  adding  binder 
decomposition  products  to  the  flame.  One  of  the  main  drivers  behind  this  work  was  that 
experimental  propellant  species  profiles  (such  as  CN)  seemed  to  be  more  extended  than 
respective  CN  profiles  in  HMX  monopropellant  models  [58].  One  of  the  goals  of  this  work  is 
to  determine  if  this  could  be  due  to  the  presence  of  the  binder.  Bum  rate  prediction  was  not  a 
goal  of  this  work.  In  this  model,  the  gases  from  the  HMX  and  the  binder  are  assumed  to  be 
premixed  instantaneously,  although  the  propellant  is  heterogeneous  (200  and  20  |im  HMX 
particles  in  the  binder).  Since  the  stoichiometries  of  the  binder/plasticizer  and  the  HMX  are 
similar,  important  diffusion  flames  are  not  formed  as  in  AP  propellants  although  there  will  be 
some  diffusional  mixing.  Diffusion  flame  modeling  incorporating  detailed  kinetics  is  being 
performed  [65],  but  is  outside  the  scope  of  this  work. 
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The  first  results  described  here  involve  the  reaction  scheme  of  Hatch  [57]  for  HMX 
combustion  which  includes  77  gas  phase  reactions  involving  HCN,  CN,  N,  H,  NH3,  NH2, 
N2,  NH,  HNO,  HNCO,  NCO,  NO2,  N20,  NO,  CH20,  HCO,  CO,  CO2,  O,  02,  OH,  H20, 
H2,  H202,  N2H,  and  H02.  This  can  be  compared  to  Melius'  131  reaction  set  which  adds 
reactions  involving  HMX,  HMX  radical  (HMX  minus  N02),  H2CNN02,  C2N2,  H2CN, 
HONO,  and  HOCN.  The  reaction  H2CN  +  M  — »  H  +  HCN  +  M  from  Melius  was  added  to 
Hatch's  reaction  set  The  second  reaction  mechanism  used  is  the  Miller  Bowman  mechanism 
for  nitrogen  chemistry  in  flames  [30],  as  extended  by  Zabamick  [11].  This  mechanism 
includes  252  reactions  involving  54  species.  In  addition  to  the  species  included  in  the  Hatch 
mechanism,  the  mechanism  includes  hydrocarbon  species  as  well  as  HCNO,  HOCN,  HONO, 
HNO3,  and  C2N2. 

The  HMX  decomposition  was  modeled  as  occurring  by  one  of  two  pathways: 

HMX  (C4H8N8O8)  ->  2  CH2O  +  2  H2CN  +  2  N2O  +  2  NO2 

HMX  -»  4  CH2O  +  2  N2  +  4/3  N2O  +  4/3  NO2 

The  first  pathway  assumes  equality  between  the  gas  and  liquid  phase  decomposition  reactions 
(described  above).  The  second  pathway  is  from  Kubota  [40].  Relative  equality  between  the 
H2CN/NO2  (gas  phase  decomposition)  and  CH2O/N2O  (liquid  phase  decomposition) 
pathways  has  been  found  experimentally  [47]  and  predicted  theoretically  [58]  for  HMX.  The 
"branching  ratio"  between  the  two  pathways  is  probably  dependent  on  pressure  and  heating 
rate.  RDX  is  found  to  break  down  mostly  through  the  H2CN/NO2  pathway  because  of  its 
higher  vapor  pressure  and  consequently  greater  tendency  to  vaporize  before  decomposition 
[47,58].  The  TMETN  in  propellant  HMXl  was  assumed  to  immediately  break  up  via: 

TMETN  (C5H9N2C>9)  — » 3  CH20  +  3  NO2  (Hatch  mechanism) 

TMETN  -*  3  CH20  +  3  NO2  +  C2H4 

Note  that  neither  scheme  is  stoichiometrically  balanced.  These  two  reactions  appear  to 
be  reasonable  assumptions  [39,57,62],  although  to  be  completely  correct  stoichiometrically, 
the  reaction  should  show  some  NO/CO/hydrocarbon/aldehyde  formation  (seen  experimentally 
in  Reference  62),  e.g.:  TMETN  (C5H9N3O9)  -*  2  CH20  +  2  N02  +  NO  +  CH2CHO.  Two 
decomposition  schemes  were  used  for  the  binder  (R-18,  -[(-CH2-CH2-0-(C=0)-(CH2)4- 
(C=0)-0-]n): 

binder  (C8H12O4)  ->  3  C2H4  +  2  C02. 

binder  — »  4  CH20  (Hatch  mechanism) 

Again,  one  of  the  reactions  is  not  stoichiometrically  balanced  because  of  the  lack  of 
hydrocarbons  in  the  Hatch  mechanism.  The  flame  combustion  mechanisms  of  species  such  as 
C2H4  are  relatively  well  understood,  at  least  in  flames  with  02  as  the  oxidizer  [63].  The 
formation  (and  incomplete  combustion)  of  solid  carbon  (soot)  may  also  be  important  [681  and 
is  not  included.  In  total  then,  there  are  four  sets  of  modeling  results,  the  two  HMX  paths  with 
the  Hatch  mechanism,  and  the  two  HMX  paths  with  the  Miller-Bowman  mechanism.  The 
major  species  results  for  these  four  sets  are  shown  in  Figures  93-96. 
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Figure  93.  Results  of  model  for  propellant  HMX1,  15  atm.  Miller-Bowman 
mechanism;  HMX  — >  4  CH2O  +  2  N2  +  4/3  N2O  +  4/3  NO2. 
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Figure  94.  Results  of  model  for  propellant  HMX1,  15  atm.  Miller- Bowman 
mechanism;  HMX  — >  2  CH2O  +  2  H2CN  +  2  N2O  +  2  NC>2* 
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To  assess  the  utility  of  the  model  predictions,  they  must  be  compared  to  experimental 
data.  For  propellant  HMX1,  the  absorption  data  of  Vanderhoff  [12],  the  CARS  data  of 
Stufflebeam  [13],  and  the  CN  LIF  profile  [4,5]  can  be  compared  with  the  model  predictions. 
As  mentioned  previously,  the  temperature  profile  used  in  the  modeling  comes  from  the 
absorption  data  of  Vanderhoff  [12].  In  CARS  experiments  at  23  atm,  Stufflebeam  measured 
the  concentration  of  N2,  CO,  H2  and  temperature  in  an  HMX1  propellant  flame  at  distances 
about  1  mm  and  greater  above  the  surface  [13].  In  absorption  experiments  at  20  atm, 
Vanderhoff  measured  the  concentration  profiles  of  CN,  NH,  and  OH  [12].  Compared  to  LIF 
experiments,  the  absorption  experiments  have  the  advantage  of  yielding  absolute 
concentrations,  rather  than  relative  concentration  profiles.  The  spatial  and  temporal  resolution 
and  sensitivity  of  LIF  is  potentially  greater  than  that  of  absorption,  however.  The  experimental 
work  at  the  Astronautics  Laboratory  has  focussed  on  measurements  of  species  and  temperature 
profiles  using  LIF.  Laser-Raman  scattering  with  a  308  nm  laser  was  attempted  for  major 
species  and  temperature  measurements,  but  was  unsuccessful  due  to  large  LIF  interferences. 
Previous  work  has  demonstrated  fairly  successful  CN  and  OH  profile  measurements  in  several 
propellants  at  pressures  up  to  3.5  MPa  (500  psi)  [4,5].  OH  LIF  temperature  measurements 
were  generally  unsuccessful  because  of  the  conflicting  requirements  of  wide  detection 
bandpass  and  high  scattered  light  rejection  (see  earlier  section),  although  wide  bandpass  OH 
LIF  temperature  measurements  were  successfully  made  at  NWC  in  atmospheric  pressure  laser- 
supported  propellant  flames  even  in  metallized  propellants  [10]. -One  common  aspect  to  all 
these  measurements  is  that  successful  data  collection  below  about  0.2  mm  above  the  surface  is 
rare,  due  to  surface  attenuation  of  laser  and  absorption  beams. 

The  published  CARS  measurements  are  of  major  species  above  about  1  mm.  These 
measurements  should  agree  with  the  model  since  they  reflect  the  equilibrium  flame  conditions. 
Any  differences  are  due  to  problems  with  the  stoichiometry  when  using  the  Hatch  mechanism 
[7,8].  However,  the  radical  measurements  should  be  sensitive  to  the  model  assumptions.  For 
example,  the  CN  measurements  are  compared  with  the  four  modeling  cases  in  Figures  97  and 
98.  In  Figure  97,  it  can  be  seen  that  the  Kubota  HMX  stoichiometry  appears  to  be  totally 
inconsistent  with  the  experimental  data  [7,8],  However,  in  Figure  98,  where  the  binder  is 
allowed  to  form  C2H4,  it  can  be  seen  that  the  Kubota  stoichiometry  actually  gives  closer 
agreement  to  the  experimental  data  than  the  other  case.  This  is  because  of  the  interaction  of  the 
CHX  species  with  the  NOx  species  in  the  Miller-Bowman  mechanism  to  form  CN  (Appendix 
2).  This  demonstrates  that  CN  formation  requires  either  the  presence  of  CHX  species  (not 
found  in  a  flame  where  the  only  fuel  is  CH2O  [52])  or  the  presence  of  HCN.  From  the 
measurements  of  Parr  [10],  it  is  apparent  that  CN  is  present  in  pure  HMX  flames,  so  it  still 
seems  unlikely  that  Kubota's  stoichiometry  is  correct.  However,  the  presence  of  CN  in 
CH2O/NO2  flames  has  led  to  speculations  about  CN  formation  through  the  reaction  of  HCO 
with  NO  and  N2O  [52],  so  this  is  still  an  area  of  uncertainty.  Another  interesting  point  that  is 
seen  in  Figures  93-96  is  that  the  Hatch  mechanism  consistently  predicts  slower  conversion  of 
species  to  final  products. 

Comparisons  of  the  NH  absorption  data  with  the  model  are  found  in  Figures  99  and 
100.  The  NH  data  is  less  informative  than  the  CN  data  because  the  predictions  of  the  various 
models  are  similar.  OH  comparisons  are  shown  in  Figures  101  and  102,  and  are  also  less 
interesting  than  the  CN  data  because  OH  acts  as  an  equilibrium  (product)  species,  rather  than  a 
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"reaction  zone"  species  as  do  NH  and  CN.  The  behavior  of  other  species  in  this  flame  is  quite 
interesting,  as  shown  in  Figures  103  (Hatch  mechanism)  and  104  (Miller-Bowman 
mechanism). 
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Figure  97.  Comparison  of  experimental  data  and  model  results  for  CN  in 
HMX1, 15  atm.  Hatch  mechanism;  "CN"=HMX  -4  2  CH20  +  2  H2CN  +  2  N20  +  2  NO2, 
"CN  K"=HMX  ->  4  CH20  +  2  N2  +  4/3  N20  +  4/3  NQ2. 
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Figure  98.  Comparison  of  experimental  data  and  model  results  for  CN  in 
HMX1, 15  atm.  Miller-Bowman  mechanism;  "CN"=HMX  -*  2  CH20  +  2  H2CN  +  2  N20 
+  2  NO*  "CN  K"=HMX  ->  4  CH20  +  2  N2  +  4/3  N20  +  4/3  NO* 
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Figure  99.  Comparison  of  experimental  data  and  model  results  for  NH  in 
HMX1, 15  atm.  Hatch  mechanism;  "NH"=HMX  -»  2  CH20  +  2  H2CN  +  2  N20  +  2  N02, 
"NH  K"=HMX  -4  4  CH20  +  2  N2  +  4/3  N20  +  4/3  N02. 
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Figure  100.  Comparison  of  experimental  data  and  model  results  for  NH  in 
HMX1, 15  atm.  Miller-Bowman  mechanism;  "NH"=HMX  ->  2  CH20  +  2  H2CN  +  2  N20 
+  2  NO*  "NH  K"=HMX  ->  4  CH20  +  2  N2  +  4/3  N2Q  +  4/3  NO* 
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Figure  101.  Comparison  of  experimental  data  and  model  results  for  OH  in 
HMX1, 15  atm.  Hatch  mechanism;  "OH"=HMX  ->  2  CH20  +  2  H2CN  +  2  N20  +  2  NO* 
"OH  K"=HMX  4  CH20  +  2  N2  +  4/3  N20  +  4/3  N02. 
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Figure  102.  Comparison  of  experimental  data  and  model  results  for  OH  in 
HMX1, 15  atm.  Miller-Bowman  mechanism;  "OH"=HMX  — »  2  CH2O  +  2  H2CN  +  2  N2O 
+  2  NO2.  "OH  K"=HMX  4  CH20  +  2  N2  +  4/3  N20  +  4/3  NO2. 
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Figure  103.  Radical  profiles  for  HMX1,  15  atm.  Hatch  mechanism;  HMX  ->  2 
CH20  +  2  H2CN  +  2  N20  +  2  NO2. 
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Figure  104.  Radical  profiles  for  HMX1,  15  atm.  Miller-Bowman  mechanism;  HMX 
2  CH20  +  2  H2CN  +  2  N20  +  2  NC^. 


HMX2  PROPELLANT 

For  propellant  HMX2,  the  probe  sampling  measurements  of  Kubota  [40]  and  the  CN 
profile  (using  LIF)  of  Edwards  [4,5]  can  be  compared  to  model  predictions.  Kubota’s  probe 
sampling  measurements  were  taken  in  two  zones  of  the  HMX2  flame  at  20  atm,  the 
"preparation  zone"  and  the  luminous  flame  [40].  Products  were  apparently  analyzed  by  GC. 
This  type  of  analysis  precludes  the  measurement  of  reactive  species  such  as  NO2,  showing 
NO,  N2,  N2O,  CO,  CO2.  H2,  and  "hydrocarbons."  It  is  not  clear  that  HCN  would  be  seen 
by  such  a  detection  system. 

The  results  for  the  four  cases  for  propellant  HMX2  are  illustrated  in  Figures  105-108. 
The  model  temperature  profile  is  taken  from  Kubota  [40].  KubQta  defined  the  relatively 
constant  temperature  region  from  0.5  to  1  mm  as  the  preparation  zone,  with  the  luminous  flame 
beginning  at  the  rapid  temperature  increase  at  about  1  mm  above  the  surface.  The  complexity 
of  the  propellant  chemistry  is  apparent,  with  several  "reaction  zones"  definable.  Near  the 
surface,  H2CN  (not  shown,  profile  indistinguishable  from  NO2  profile)  and  NO2  rapidly 
react,  forming  HCN  and  NO.  CH2O  reacts  less  rapidly,  persisting  out  to  0.5  mm  with  the 
Hatch  mechanism.  The  most  striking  feature  is  the  lack  of  HCN  reaction  predicted  in  Figures 
105  and  107  for  the  HMX  — »  2  CH20  +  2  H2CN  +  2  N20  +  2  N02  pathway,  even  at  the 
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adiabatic  flame  temperature.  Under  these  flame  conditions,  HCN  is  less  reactive  than  NO, 
which  is  consumed  in  the  luminous  flame.  It  has  generally  been  found  that  HCN  is  less 
reactive  than  CH2O  in  NO2  oxidized  environments  [66].  If  correct,  this  large  (in  both  spatial 
and  absolute  concentration)  HCN  level  should  be  accessible  with  CARS  [13]. 

There  are  two  data  sets  for  propellant  HMX2  that  can  be  used  for  model  testing.  The 
first  set  is  the  probe  results  of  Kubota  [40];  the  second  data  set  is  the  CN  LEF  profile  [4,5].  In 
Kubota's  experiment,  GC  (?)  samples  were  extracted  from  the  luminous  flame  and  the 
"preparation  zone"  (the  nonluminous  zone  between  the  surface  and  the  luminous  flame)  of  an 
HMX  flame  at  20  atm.  The  two  zones  can  be  seen  in  the  temperature  profile  as  the  flat  regions 
at  about  0.5  mm  above  the  surface  (preparation  zone)  and  above  1.5  mm  above  the  surface 
(luminous  flame)  at  30  atm.  The  probe  data  are  at  20  atm  and  the  temperature  profile  was 
measured  at  30  atm.  For  comparison,  I  assumed  probe  measurements  at  30  atm  would  be  the 
same  as  those  at  20  atm.  Kubota's  results  are  shown  in  Table  13.  The  luminous  flame  results 
are  very  close  to  the  predicted  equilibrium  species  concentrations,  as  is  the  temperature 
measured  by  microthermocouples.  In  the  preparation  zone,  Kubota  found  large  amounts  of 
NO  and  N2O,  as  well  as  "hydrocarbon  fragments  and  solid  carbon  generated  at  the  burning 
surface"  [40].  Since  the  Hatch  mechanism  described  in  this  paper  doesn’t  include 
hydrocarbons  (CH2O  is  not  a  "hydrocarbon"  in  the  strict  sense  of  the  word),  matching  this  last 
observation  is  impossible  for  this  mechanism.  However,  it  was  also  found  that  the  model 
could  not  match  several  of  the  other  observations  for  either  mechanism.  The  only  adjustable 
parameters  in  the  model  are  the  starting  species.  Two  inconsistencies  were  found.  First,  when 
significant  amounts  of  HCN  were  formed  by  assuming  the  H2CN/CHO2/NO2/N2O  pathway  in 
the  model,  the  HCN  persisted  even  through  the  luminous  flame.  Second,  if  the 
CH2O/N2O/NO2/N2  (Kubota)  pathway  was  assumed,  then  very  little  NO  was  predicted  in  the 
preparation  zone  since  the  CH2O/N2O/NO2/N2  pathway  doesn’t  form  large  amounts  of  NO 
compared  to  the  other  pathway.  Several  explanations  for  these  discrepancies  are  possible. 
First,  the  slow  reaction  of  HCN  at  this  propellant's  relatively  low  flame  temperature  may  be 
due  to  missing  HCN-consuming  reactions  in  the  mechanism  (although  both  mechanisms  show 
the  lack  of  HCN  reaction).  At  the  higher  temperatures  found  in  propellant  HMX1,  HCN  is  not 
predicted  to  survive  the  luminous  flame  (as  shown  in  the  previous  section).  Second,  Kubota's 
temperature  profile  may  be  wrong.  Microthermocouple  measurements  in  propellant  flames  are 
not  trivial,  although  it  seems  likely  that  die  the  observation  of  two  zones  in  these  slow-burning, 
fuel-rich  propellants  is  correct  TTiird,  it  is  conceivable  that  the  gas  sampling  system  may  have 
destroyed  HCN  and  NO2  present  in  the  sample,  thus  incorrectly  implying  their  absence. 
Fourth,  the  initial  species  assumptions  may  be  wrong. 

One  problem  with  the  Kubota  stoichiometry  is  that  the  mechanism  forming  the  products 
is  not  clear.  The  H2CN/NO2  mechanism  (gas  phase  HMX  decomposition)  is  supported  by 
quantum  mechanical  calculations  [59]  and  experiments  [47],  while  the  condensed  phase 
CH2O/N2O  pathway  is  supported  by  many  slow  heating  rate,  condensed  phase  experiments 
[67]  and  the  presence  of  CH2O  and  N2O  in  high  thermolysis  experiments  [47].  HCN  has  been 
detected  in  HMX  propellant  combustion  with  CARS,  although  the  amount  has  not  yet  been 
quantified  [13].  TTius  it  seems  unlikely  that  CH2O  is  the  only  fuel  molecule  formed.  It  also 
seems  likely  that  the  high  pressures  of  the  propellant  combustion  involve  extrapolating  the 
known  rate  constants  well  beyond  their  validated  range.  Thus  the  lack  of  agreement  between 
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the  experimental  data  and  the  model  is  not  necessarily  an  indication  that  the  model  is  wrong, 
but  rather  indicates  where  improvements  in  understanding  the  model  and  experiments  are 
needed  One  partial  confirmation  of  the  model  is  found  in  CN  LIF  profiles  in  propellant 
HMX2.  As  shown  in  Figures  109  and  110,  the  experimental  data  shows  a  broad  (but  weak) 
CN  peak  at  a  relatively  large  distance  from  the  propellant  surface  [4,5].  The  models  show  that 
even  though  large  amounts  of  HCN  are  present  in  the  model  flame,  the  CN  concentration  is 
very  low  because  of  the  low  temperature.  The  CN  that  is  present  is  "equilibrium"  CN  in  the 
luminous  flame,  rather  than  "reaction  zone"  CN  as  is  normally  seen  in  flames.  The 
experimental  data  appears  to  confirm  this,  both  in  the  low  intensity  (concentration)  and  lack  of 
a  peak  near  the  surface.  However,  the  data  has  a  very  low  signal-to-noise  ratio,  so  this  is  only 
a  partial  confirmation  until  better  data  is  collected 
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Figure  105.  Results  of  model  for  propellant  HMX2,  30  atm.  Hatch  mechanism; 
HMX  2  CH2Q  +  2  H2CN  +  2  N2Q  +  2  NC^  (case  #1). 
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Figure  106.  Results  of  model  for  propellant  HMX2,  30  atm.  Hatch 
HMX  -»  4  CH20  +  2  N2  +  4/3  N20  +  4/3  NO2  (case  #2) . 
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Figure  107.  Results  of  model  for  propellant  HMX2,  30  atm.  Miller-Bowman 
mechanism;  HMX  -»  2  CH20  +  2  H2CN  +  2  N20  +  2  NC^  (case  #3). 
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Figure  108.  Results  of  model  for  propellant  HMX2,  30  atm.  Miller-Bowman 
mechanism;  HMX  -»  4  CH2O  +  2  N2  +  4/3  N2O  +  4/3  NO2  (case  #4). 


Table  13.  Kubota  probe  results  for  prop.  HMX2  (20  atm)  (mole  fractions,  H20-free  basis) 

[40]. 

■PREP.  ZONE"  "LUMINOUS  FLAME" 


0.5  mm 

1.75  mm 

SPECIES 

EXPT 

EXPT 

EQ  CODE 

NO 

0.210 

— — 

— - 

n2 

0.098 

0.228 

0.245 

n2o 

0.068 

— 

— 

CO 

0.180 

0.422 

0.420 

CO2 

0.073 

0.041 

0.032 

h2 

0.067 

0.267 

0.303 

"HC  +  Cs" 

0.304 

— 

— 

T,  K 

1280 

1870 

2080 
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Table  14.  Model  results  at  30  atm  for  propellant  HMX2  (mole  fractions,  H20-free  basis). 
Case  #*s  listed  in  legends  of  Figures  105-108. 


Case  #1 

Case  #2 

Case  #3 

Case  #4 

Species 

0.5  mm 

1 .75  mm 

0.5  mm 

1 .75  mm 

Oj^nm 

iT^nm 

0.5  mm 

1 .75  mm 

NO 

0.073 

0.006 

0.087 

0.015 

0.143 

0.029 

0.086 

0.036 

n2 

0.030 

0.273 

0.171 

0.307 

0.038 

0.240 

0.166 

0.265 

N20 

0.188 

— 

0.100 

— 

0.160 

— 

0.085 

.... 

CO 

0.318 

0.271 

0.287 

0.373 

0.202 

0.267 

0.252 

0.316 

C02 

0.003 

0.149 

0.001 

0.077 

0.087 

0.122 

0.110 

0.089 

h2 

0.158 

0.141 

0.184 

0.228 

0.070 

0.141 

0.199 

0.204 

HCN 

0.167 

0.159 

— 

— 

0.175 

0.155 

0.0002 

0.019 

CH20 

0.062 

~ 

0.167 

— 

0.0001 

— 

0.0006 

.... 

HC* 

— 

— 

— 

— 

0.070 

0.028 

0.064 

0.079 

*  primarily  C2H4,  C2H2 
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Figure  109.  Comparison  of  CN  LIF  profile  and  modeling  data.  HMX2,  30  atm. 
Hatch  mechanism. 
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Figure  110.  Comparison  of  CN  LIF  profile  and  modeling  data.  HMX2,  30  atm. 
Miller- Bowman  mechanism. 


AN1  PROPELLANT 

Recently,  temperature  profiles  for  an  AN  propellant  flames  have  became  available  [69]. 
This  allows  chemical  kinetic  modeling  of  the  AN  propellant  flame  chemistry,  similar  to  that 
performed  for  HMX  propellants.  There  is  less  experimental  data  available  for  AN  propellants, 
due  to  experimental  difficulties.  AN  propellant  flames  are  very  "sooty,"  and  are  thus  difficult 
to  probe  with  laser  diagnostics  [4,5].  Thus,  comparing  these  modeling  predictions  with 
experimental  data  will  be  difficult 

The  temperature  profile  used  was  taken  from  microthermocouple  measurements  [69]. 
The  experimental  data  for  15  atm  is  shown  in  Figure  111.  The  data  in  Figure  11 1  is  from 
thermocouples  made  from  75  |im  wire  (Pt/Pt-13%Rh)  with  a  bead  size  of  approximately  200 
pm.  Measurements  with  smaller  wire  are  planned.  This  is  the  same  size  wire  as  that  used  by 
Kubota  for  his  published  temperature  profiles  in  HMX  propellants  [40].  Smaller 
thermocouples  may  indicate  a  more  rapid  temperature  rise  above  the  surface,  although  the 
surface  temperature  and  final  flame  temperature  probably  will  not  change  significantly.  The 
measured  surface  temperature  of  600  K  is  in  line  with  other  AN  burning  surface  temperature 
measurements  [74], 
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The  first  reaction  set  used  was  taken  from  a  paper  by  Hatch  [27],  which  includes  77 
reactions  involving  26  species.  An  additional  reaction  (HNO3  -*  NO2  +  OH)  was  taken  from  a 
paper  by  Zabamick  [11].  The  propellant  modeled  was  an  energetic  binder  AN  propellant 
("AN1":  67  wt%  AN  (50%  15  *tm,  50%  190  [im),  21  wt%  GAP,  12  wt%  TMETN),  with  a 
burning  rate  of  0.3  cm/s  at  15  atm  (0.5  g/cm2s).  The  following  reactions  were  assumed  to 
happen  at  the  propellant  surface  for  the  chemical  kinetic  calculations: 

AN:  NH4NO3  -*■  NH3  +  HNO3 

GAP:  C3H5N3O  ->  CH20  +  HCN  +  N2 

TMETN:  C5H9N3O9  -♦  3  CH20  +  3  NC^ 

These  reactions  are  the  "baseline"  case.  Note  that  the  Hatch  mechanism  doesn't  include 
hydrocarbons,  so  the  GAP  and  TMETN  reactions  are  not  stoichiometrically  balanced.  The 
baseline  stoichiometry  for  the  propellant  is  then  (mole  fractions): 

AN1  -+  0.323  NH3  +  0.323  HNO3  +  0.082  HCN  +  0.082  N2  +  0.136  CH20  +  0.054  NO2. 
The  baseline  reactions  were  postulated  based  mainly  on  high  rate  thermolysis  data  of  Brill 
[62,70,71].  The  model  is  a  premixed  model  so  no  diffusion  flame  effects  are  currently 
modeled.  Detailed  chemical  kinetic  modeling  of  diffusion  flames  is  still  in  its  infancy  [65]. 

Species  profiles  for  the  baseline  stoichiometry  with  the  Hatch  mechanism  are  shown  in 
Figures  112  and  113.  The  model  flame  structure  is  complex,  with  several  "reaction  zones" 
discemable.  HNO3,  CH20,  and  part  of  the  NH3  are  seen  to  react  very  close  to  the  surface, 
with  the  NO2  formed  from  the  HNO3  reacting  with  the  remaining  NH3.  HCN  was  not 
predicted  to  react  very  readily  at  the  low  temperatures  of  the  AN1  propellant  flame.  This  is 
consistent  with  the  model  predictions  in  HMX  propellant  flames.  However,  the  reaction  set 
may  not  include  important  HCN  consuming  reactions.  The  addition  of  hydrocarbon  reactions 
to  die  mechanism  will  allow  better  predictions,  especially  for  the  GAP  binder  where  Brill  found 
significant  amounts  of  hydrocarbons  formed  during  thermolysis  [70]. 

For  comparison,  the  profiles  for  a  pure  AN  flame  (NH3/HNO3M3.5/O.5)  of  the  same  T 
profile  and  mass  burning  rate  is  shown  in  Figures  1 14  and  1 15.  The  profiles  for  NH3,  HNO3, 
and  NO2  are  similar  to  those  in  the  propellant  Interestingly,  the  equilibrium  O2  concentration 
(0.14  mole  fraction)  is  not  reached,  with  relatively  equal  amounts  of  O2  and  NO  formed  as 
final  products.  At  15  atm,  the  equilibrium  conditions  are  Tfi=1245  K, 
H20/02/N2=0.57/0. 14/0.29  mole  fractions. 

As  a  rough  approximation  to  the  diffusion  flame  behavior  in  the  propellant,  one  could 
assume  that  the  fast  reactions  involving  HNO3  and  NH3  occur  in  the  monopropellant  (AN) 
flame,  with  the  "leftovers"  from  that  flame  reacting  with  the  binder-generated  species.  For 
example,  assuming  that  the  species  profiles  at  0.55  mm  in  Figure  1 14  are  the  products  from  the 
AN  monopropellant  flame  (NH4NO3  -*  1.2  H20  +  0.25  N2  +  0.6  NH3  +  0.9  NC^),  the 
chemistry  of  the  AN/binder  diffusion  flame  could  be  roughly  approximated  as  the  flame 
between  the  AN  monopropellant  flame  products  and  the  binder  pyrolysis  products.  This  leads 
to  a  flame  (approximated  as  premixed)  with  the  following  stoichiometry: 
N02/CH20/N2/H20/HCN/NH3  =  0.264/0.104/0.124/0.297/0.063/0.148.  The  species 
profiles  from  such  a  flame  are  shown  in  Figures  116  and  117.  Reactions  between  the 
monopropellant  flame  products  and  the  binder  products  do  not  begin  until  the  maximum  flame 
temperature  is  nearly  reached.  This  indicates  that  the  assumption  that  the  monopropellant 
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products  and  the  binder  products  have  time  to  mix  may  not  be  unrealistic.  This  stoichiometry 
may  be  compared  with  the  "secondary  reaction"  products  seen  by  Brill  from  AN  thermolysis 
(HNO3  and  NH3  are  the  initial  products  seen) :  NH4NO3  -*■  0.68  H2O  +  0.24  NH3  +  0.6  H2 
+  0.96  NO2  +  0.4  N2O.  This  stoichiometry  was  estimated  from  Brill’s  data  which  did  not 
include  H2O,  H2,  N2,  or  O2. 

Alternatively,  one  could  assume  that  the  AN  reacted  mostly  in  the  liquid  phase: 
NH4NO3  -*  2  H2O  +  N2O.  Then  the  gas  phase  reactions  would  be  mostly  fuel/N20  reactions. 
The  model  predictions  for  this  situation  are  shown  in  Figures  118  and  119.  The  behavior  is 
similar  to  that  shown  in  Figures  116  and  117  for  the  monopropellant  products/binder  flame, 
with  little  reaction  before  about  0.5  mm  above  the  surface  (T-1200  K).  This  is  because  the 
model  predicts  that  N2O  is  less  reactive  at  low  temperatures  than  HNO3. 


AN  Propellant,  200  psi 


Figure  111.  AN1  Propellant  Temperature  Profile  -  15  atm.  For  model,  gas  phase 
T  profile  approximated  as  linear  between  x=0,  T=600  K  and  x=0.9  mm,  T=1700  K. 
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Figure  112.  Model  results  for  baseline  case,  major  species. 


Figure  113.  Model  results  for  baseline  case,  selected  minor  species. 
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Figure  114.  Model  results  for  NH3/HNO  3=0.570.5,  major  species. 
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Figure  115.  Model  results  for  NH3/HN03=0.5/0.5,  minor  species. 
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Figure  116.  Model  results  for  flame  between  monopropellant  products  and 
binder  pyrolysis  products,  major  species. 
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Figure  117.  Model  results  for  flame  between  monopropellant  products  and 
binder  pyrolysis  products,  minor  species. 
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Some  very  recent  modeling  of  the  AN1  propellant  has  been  done  with  the  extended 
Miller-Bowman  mechanism  (described  in  the  previous  HMX  propellant  discussions).  This 
mechanism  includes  hydrocarbons,  so  the  stoichiometry  for  the  various  ingredients  is 
somewhat  different  (but  still  based  on  the  data  of  Brill): 

AN:  NH4NO3  -*■  NH3  +  HNO3 

GAP:  C3H5N3O  -►  CH4  +  HCN  +  N2  +  CO 

TMETN:  C5H9N3O9  -*  3  CH2O  +  3  NO2  +  C2H4 

The  initial  results  for  this  case  are  shown  in  Figure  120  (major  species)  and  Figure  121  (minor 
species).  The  temperature  profile  is  again  that  shown  in  Figure  111.  It  is  evident  from  these 
figures  that  the  Miller-Bowman  mechanism  predicts  much  more  complete  reaction  close  to  the 
surface  that  does  the  Hatch  mechanism.  No  intermediate  HCN  "plateau"  is  seen,  with  rapid 
reactions  occurring  when  the  flame  temperature  reaches  about  800  K.  For  comparison,  the 
species  profiles  for  a  NH3/HNO3  flame  with  the  same  temperature  profile  are  shown  in  Figures 
122  and  123.  Again,  reactions  are  rapid,  with  complete  consumption  of  fuel  and  oxidizer 
within  0.4  mm  of  the  surface.  These  rapid  reactions  cast  some  doubt  on  the  multiple  reaction 
zones  seen  with  the  Hatch  mechanism.  For  completeness,  models  with  the  Hatch  1 
stoichiometry  should  be  run  with  the  Miller-Bowman  mechanism.  At  the  present  time,  this  has 
not  been  done.  It  may  be  that  the  hydrocarbons  in  the  mechanism  are  responsible  for  the  rapid 
reactions,  rather  than  the  presence  of  more  HCN-consuming  reactions. 
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Figure  120.  Model  results  (major  species)  for  AN1  propellant  flame,  Miller- 
Bowman  mechanism. 
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Figure  121.  Model  results  (minor  species)  for  AN1  propellant  flame,  Miller- 
Bowman  mechanism. 
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Figure  122.  Model  results  (major  species)  for  AN  in  AN1  propellant  flame, 
assuming  AN  -+  HNO3  +  NH3. 


Figure  123.  Model  results  (minor  species)  for  AN  in  AN1  propellant  flame. 
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Modeling  of  AN  propellant  flames  using  "global"  kinetics  has  been  performed  by 
Beckstead  [74]  using  the  Separate  Surface  Temperature  (SST)  model,  the  most  up-to-date 
version  of  the  Beckstead-Derr-Price  (BDP)  model.  AN  will  not  bum  as  a  monopropellant  at  15 
atm  so  the  monopropellant  flame  is  probably  not  an  important  pan  of  the  propellant  combustion 
chemistry.  At  15  atm,  the  "monopropellant  flame  height"  can  be  found  by  extrapolation  to  be 
about  700  pm  [74].  Thus  the  combustion  is  controlled  by  the  diffusion  flame  between  the  AN 
and  the  binder  (the  "primary"  flame).  Beckstead  calculates  the  diffusional  distance  for  the 
primary  flame  and  multiplies  it  by  3  to  obtain  the  "primary  flame  height."  The  diffusional 
distance  is  calculated  by  a  numerical  approximation  to  a  Burke  Schumann  calculation: 
x*/bi=(yst/al)a2(afi/aoxi)a3 

where: 

bi=distance  from  center  of  oxidizer  particle  i  to  center  of  binder  layer  surrounding 

particle 

ySt=stoichiometric  oxidizer  to  fuel  ratio  (mass) 
ag=mass  fraction  of  binder  associated  with  oxidizer  fraction  i 
otoxi=ntass  fraction  of  oxidizer  associated  with  oxidizer  fraction  i 
al,a2,a3=0.834, 0.858,1. 871  (flame  closes  over  binder) 
al,a2,a3=0.472, 0.929, 1.604  (flame  closes  over  oxidizer) 

For  AN  propellants,  the  stoichiometric  ratio  was  calculated  assuming  the  final  diffusion  flame 
products  were  CO  and  H2O.  Defining  the  Under  as  TMETN  (C5H9N3O9)  +  GAP  (C3H5N3O) 
yields  a  fuel  C0.87H1.48N0.7gO0.64-  Assuming  the  AN  contributes  0.5  mole  02/mole  AN 
(NH4NO3  ->  N2  +  2  H2O  +  0.5  62),  the  stoichiometric  ratio  for  AN  is  0.425  gm  binder/gm 
AN  (1.03  mole  binder/mole  AN),  or  ySt=2.35  gm  AN/gm  binder.  In  other  words,  NH4NO3 
+  1.03  Cq.87Hi.48No.780o.64  0.90  CO  +  1.40  N2  +  2.76  H2O.  yst  is  about  10  for  non- 

oxygenated  (HTPB)  binders  [74].  The  GAP/TMETN  binder  could  also  be  modeled  as  an 
energetic  binder  with  a  separate  flame.  However,  heat  of  explosion  (Hex)  calculations,  with  an 
assumed  stoichiometry  of  binder  -+  0.39  N2  +  0.64  H2O  +  0.11  H2  +  0.87  Cs  (solid  carbon), 
yield  a  Hex  value  of  -800  cal/g.  Assuming  the  products  are  CO,  H2  and  Q  yields  a  HeX  of 
-550  cal/g.  The  800  cal/g  result  fits  a  correlation  previously  seen  between  Hex  and  the 
adiabatic  flame  temperature  [5],  so  the  800  cal/g  value  is  probably  more  correct.  With  this  low 
Hex  value,  it  was  assumed  that  the  energetic  binder  reactions  between  the  HCN,  CH2O  and 
NO2  were  slow  enough  that  mixing  between  the  AN  and  the  binder  products  occurred  before 
the  energetic  binder  reactions.  Chemical  kinetic  modeling  of  possible  energetic  binder  reactions 
could  not  be  made  because  of  the  lack  of  reactions  involving  Cs  in  the  mechanism.  It  is 
possible  that  reactions  involving  Cs  are  important  in  AN  propellant  flames,  so  a  better  kinetic 
model  should  attempt  to  include  reactions  involving  Cs.  afi  and  ctoxi  were  calculated  by 
assuming  that  the  small  AN  particles  bum  with  a  stoichiometric  flame,  with  the  large  particles 
burning  the  remaining  fuel.  This  leads  to  (ccfi/ctoxi)  15^=0-425,  (afi/aoxi)i90um=0-56.  The 
amount  of  binder  associated  with  each  particle  is  needed  to  calculate  b: 
b=Di’/2[  1  -KpfiApoxi]0*5. 

where  Di'^Dj/S)0-5,  a  statistical  average  of  the  particle  diameter  (Di)  which  changes  as  the 
particle  bums,  and  <pfi  and  (poxi  are  the  volume  fractions  of  the  binder  and  AN,  respectively, 
calculated  from  the  densities  of  binder  and  AN  (1.4  and  1.7  g/cm3)  and  the  propellant 
composition.  These  calculations  lead  to  bi5^m-9.3  pm  and  bi9o^m=130  pm.  The  diffusional 
distances  are  then  calculated  assuming  the  flame  closes  over  the  oxidizer  for  the  15  pm  particles 
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and  over  the  binder  for  the  190  |im  particles  [74].  The  results  are  x*i5^m=10  pm  and 
x*190pra=l  10  pm.  In  the  SST  model,  these  diffusional  distances  are  multiplied  by  3  to  obtain 
the  "flame  height".  In  Figures  112  and  113,  the  premixed  "flame  height"  can  be  estimated  as 
300-400  pm  (very  dependent  upon  the  T  profile  used  [7]),  whereas  the  Figures  120  and  121 
yield  a  flame  height  of  about  200  pm.  From  these  calculations,  it  appears  unreasonable  to 
approximate  the  AN  propellant  flame  as  a  premixed  AN/binder  flame  at  15  atm  since  the 
diffusional  distance  is  of  the  same  order  of  magnitude  as  the  premixed  "flame  height”,  i.e.,  the 
AN  and  binder  species  must  diffuse  together  before  reacting.  This  analysis  could  also  be  made 
for  the  more  reactive  AP  propellant  flames.  For  example,  SST  calculations  for  a  stoichiometric 
AP/HTPB  diffusion  flame  yield  a  diffusional  distance  of  x*/bi=0.14  for  large  AP  particles 
(flame  closes  over  binder,  Yst=5  (CO),  <Xfi/Ooxi=0-15),  or  about  15  pm  (maximum)  for  200 
pm  AP  particles.  Chemical  kinetic  calculations  involving  fuel/QOx  reactions  would  have  to  be 
made  to  determine  if  the  kinetics  or  the  diffusion  controls  AP/HTPB  primary  flames  at  low 
(»15  atm)  pressures.  Note  also  that  some  liquid  phase  mixing  may  occur  in  low-melting, 
slow-burning  AN  propellants  [74].  Several  factors  may  contribute  to  die  much  smaller  particle 
size  influence  on  bum  rate  for  AN  propellants  versus  AP  propellants. 

If  future  temperature  measurements  with  smaller  thermocouples  determine  that  the 
temperature  profile  of  Figure  1 1 1  is  incorrect,  with  the  correct  temperature  profile  reaching  the 
maximum  closer  to  the  surface,  then  the  premixed  flame  approximation  becomes  more 
unreasonable.  Model  predictions  using  temperature  profiles  that  are  "compressed”  closer  to  the 
propellant  surface  show  that  the  species  profiles  are  very  similar  to  Figures  1 12-1 19,  except 
that  the  profiles  are  also  compressed  in  the  same  ratio  as  the  temperature  profile.  In  other 
words,  the  "reaction  zone"  is  still  close  to  the  point  where  the  temperature  profile  reaches  its 
maximum.  This  is  similar  to  the  behavior  found  in  HMX  propellant  modeling. 

There  is  not  a  lot  of  experimental  data  available  on  species  profiles/concentrations  in 
AN  flames  or  AN  propellant  flames,  especially  under  combustion  conditions.  In  addition  to 
Brill's  high  rate  thermolysis  data  [71],  there  is  some  laser-pyrolysis/gas  analysis  data  [72].  In 
the  laser  pyrolysis  experiments,  gas  analysis  from  pyrolyzed  samples  (200  cal/cm2-s,  50  psi 
He)  of  AN  yielded  roughly  equal  amounts  of  NO  and  N2,  with  relatively  little  N2O.  In 
contrast,  most  earlier  work  found  large  quantities  of  N2O  as  a  product  [72].  Apparently,  the 
products  of  AN  decomposition  are  dependent  upon  pressure/heating  rate,  with  high  heating 
rates/low  pressures  yielding  NH3/HNO3  and  low  rates/high  pressures  yielding  N2O/H2O. 
Interestingly,  the  industrial  production  of  N2O  involves  the  (slow)  heating  of  pure  AN  to  200 
C  [73].  Thus,  the  composition  of  the  gases  just  above  the  surface  of  AN  propellants  is  still  in 
doubt  (as  it  is  for  HMX  propellants). 

These  model  predictions  can  be  used  to  assess  the  chances  for  success  of  species 
measurements  or  to  suggest  good  diagnostic  targets.  Successful  optical  diagnostic 
measurements  have  been  made  in  high  pressure  HMX  propellant  flames,  using  techniques  such 
as  LIF  [4,5,7],  absorption  [12],  and  CARS  [13].  The  radical  measurements  are  difficult,  with 
relatively  low  signal-to-noise  ratios  in  the  HMX  propellant,  and  little  success  has  been  reported 
in  the  more  opaque  AN  propellant  flames  [4,5].  In  Table  15,  the  model  predictions  of  the 
concentration  of  diagnostically  accessible  (with  absorption  or  LIF)  species  in  AN  propellant 
flames  for  various  assumed  stoichiometries  are  compared  with  the  model  predictions  for  the 
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HMX  propellant  flame  at  15  atm  (for  the  Hatch  mechanism).  If  the  concentration  of  the 
particular  species  is  much  less  than  in  the  HMX  propellant  flame,  the  diagnostic  results  will  be 
difficult  to  obtain.  From  Table  15,  it  can  be  seen  that  there  is  no  easy  target  for  species 
measurements  in  the  AN  propellant  flame.  The  radicals  (OH,  CN,  NH,  and  NO)  have 
concentrations  in  the  AN  propellant  flame  at  best  equal  to  that  in  the  HMX  propellant  flame. 
This  is  due  to  the  low  temperature  and  fuel-rich  nature  of  AN  propellant  combustion.  Not 
coincidentally,  these  same  two  characteristics  of  AN  propellant  flames  lead  to  poor  combustion 
performance  (such  as  combustion  efficiency  in  aluminized  AN  propellants)  [74].  One 
interesting  difference  is  in  the  O2  concentration,  however.  O2  is  accessible  by  LIF  [75],  and 
measurements  of  O2  concentration  in  AN  propellant  flames  might  help  differentiate  between  the 
various  possible  AN  decomposition  pathways.  CARS  measurements  of  major  species  such  as 
N2O  and  NO  would  also  be  useful. 


Table  15.  Comparison  of  model  predictions  for  HMX  and  AN  propellants  at  15  atm.*  Hatch 

mechanism. 


HMX  [2] 

■  mSONHM 

AN  products/ 
binder  flame 

[CN]max,  #/cm3 

■KETCH 

2X1014 

1  X  1015 

HKETCRH 

OH 

2  X  1016 

5  X  1016 

■f«  1 

NH 

2X1015 

2X1015 

2  X  IOI4 

j  SB 

NO 

2X1019 

5  X  1018 

3  X  1018 

8  X  1018 

_ O2 _ 

9  X  1014 

2  X  1018 

5  X  1016 

1  X  1018 

*  Concentrations  shown  are  the  maximum  concentrations  predicted  by  the  model.  The  HMX 
propellant  is  73%  HMX,  17%  TMETN,  and  10%  polyester  binder  ("HMX1")  [4,5,7,12,13]. 
The  "AN  products/binder  flame"  is  discussed  for  Figures  1 14  and  1 15.  The  conversion  from 
mole  fraction  of  "A"  to  absolute  concentration  (molecules/cm3,  or  #/cm3)  is  done  for  each 
model  point  using  the  computed  density  (temperature  dependent):  mole  friction  (moles  A/moles 
gas)  X  6.023X1023  (molecules  A/moles  A)  X  gas  density  (gm  gas/cm3)  /  MWavg  (gm 
gas/mole  gas)  =  #/cm3.  For  AN  propellants,  MWavg  was  estimated  as  23,  an  average  of  H2O 
and  N2/CO. 
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SUMMARY  AND  CONCLUSIONS 

Species  data  in  propellant  flames  is  difficult  to  obtain,  especially  at  higher  pressure. 
However,  such  data  is  slowly  becoming  available,  as  shown  in  Table  16.  Of  the  reactive 
(radical)  species,  only  the  diatomics  CN,  NH,  OH,  and  NO  and  the  triatomic  NO2  (tentative) 
have  been  measured.  Other  important  molecules,  such  as  HNO  and  HCO,  have  been  measured 
under  specialized  conditions  but  are  probably  not  measurable  in  propellant  flames.  Stable 
species  in  high  concentration  that  have  been  measured  include  H2,  CO,  and  N2,  with  HCN 
identified.  Thus,  some  data  is  available.  This  data  can  be  used  to  validate  or  improve  chemical 
kinetic  models  of  propellant  flames.  There  is  also  data  available  under  noncombustion 
conditions  that  is  very  relevant  [47,62].  One  problem  area  is  that  kinetic  mechanisms  relevant 
to  propellant  combustion  have  not  been  tested  in  appropriate  flames  to  determine  the  accuracy 
of  the  parts  before  inclusion  in  the  complete  mechanism.  Thus,  when  experiment  differs  from 
model,  the  problem  can  either  be  in  the  assumptions  of  the  model  (e.g.,  initial  gas  phase 
species)  or  in  the  reaction  mechanism  itself  (e.g.,  missing  NO2  reactions).  Thus  it  is  important 
to  test  the  mechanism  in  well-controlled  flames  before  inclusion  in  the  propellant  model.  Some 
mechanism  testing  of  this  type  is  underway  (Table  17),  but  more  needs  to  be  done.  It  may  well 
be  that  many  reaction  rates  need  to  (re)measured  before  successful  propellant  flame  chemical 
kinetic  models  are  created. 


Table  16.  Propellant  flame  species  measurements.  Most  successful  measurements  usually  in 

nitramine  or  nitramine  composite  propellants. 


Author 

Diagnostic 

Reference 

Parr,  Han  son-Parr 

PUF 

CN,  NH,  OH,  NO, 
N02,T 

10 

Stuffiebeam 

CARS 

n2,  CO,  h2,  t 

13 

Vanderhoff 

Absorption 

CN,  NH,  OH,  T 

12 

Edwards 

LIF 

CN,  NH,  OH,  T 

4 

Kubota 

Probe 

NO,  N2,  N20,  CO, 
C02,  H? 

40 

Table  17.  Propellant-related  flame  measurements/modeling 


Author 

Flame 

Reference 

Zabamick 

CH4/NO2/O2 

11 

Thome 

HCN/NO2,  C2N2/NO2 

53 

Branch 

CH4/NO2/O2.  CH2O/NO2, 
CH?0/N70,  CH4/N7O 

52 
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APPENDIX  1.  HATCH  MECHANISM  [57] 
Rate  constants  for  reactions  listed  in  form  k=AT*Bexp(-E/RT) 


Reaction 

A 

R 

£ 

NH3+M-NH2+H+M 

.  920E16 

0. 

84800 

HNO+M-H+NO+M 

. 300E17 

0. 

48680 

N2H+M-N2+H+M 

.  200E15 

0. 

20000 

H2+M-H+H+M 

. 223E13 

0.5 

92600 

02+M-O+O+M 

.  185E12 

0.5 

95560 

N20+M-N2+0+M 

.  162E15 

0. 

51600 

CH20+M-HC0+H+M 

.  331E17 

0. 

81000 

HCO+M-H+CO+M 

.  145E15 

0. 

19000 

N02+M-N0+0+M 

.  110E17 

0. 

66000 

H202+M-0H+0H+M 

.  120E18 

0. 

45500 

NCO+M-N+CO+M 

.310E17 

-0.5 

48000 

NCO+OH-NO+CO+H 

, 100E14 

0. 

0 

NH3+H»NH2+H2 

.  246E14 

0. 

17071 

NH3+0-NH2+0H 

. 150E13 

0. 

6040 

NH2+OH-NH+H20 

.  125E14 

0. 

2200 

NH3+0H-NH2+H20 

.  325E13 

0. 

2120 

H+HNO-NH+OH 

.  200E12 

0.5 

1300 

HNO+OH-NO+H20 

.  360E14 

0. 

0 

NH2+HNO-NH3+NO 

. 500E14 

0. 

1000 

NH2+N0-N2H+0H 

.  468E20 

-2.46 

1876 

NH2+N0-N2+H20 

.  702E20 

-2.46 

1876 

NH+02-HN0+0 

.  112E12 

0. 

3250. 

N2H+NO-N2+HNO 

. 500E14 

0. 

0. 

N2H+OH-N2+H20 

. 300E14 

0. 

0. 

NH2+NH2-NH3+NH 

. 630E13 

0. 

10000. 

CO+OH-C02+H 

. 151E8 

1.3 

-758. 

H2+0H-H20+H 

.  520E14 

0. 

6500. 

H+02-OH+O 

.  719E17 

-0.861 

16523. 

0+H2-0H+H 

. 180E11 

1. 

8826. 

0H+0H-0+H20 

.  170E7 

2.03 

-1190. 

N20+0-N0+N0 

.  100E15 

0. 

28200. 

N20+O-N2+O2 

.  100E15 

0. 

28200. 

N20+H-N2+0H 

.760E14 

0. 

15200. 

N20+NH-N2+HN0 

.  100E12 

0.5 

3000. 

CH20+0H-HC0+H20 

. 7S3E13 

0. 

167. 

CH20+H-HC0+H2 

. 331E15 

0. 

10500. 

CH20+0-HCO+OH 

. 501E14 

0. 

4600. 

HC0+0H-C0+H20 

.  100E15 

0. 

0. 

HCO+H-CO+H2 

. 200E15 

0. 

0. 

HCO+O-OH+CO 

.  100E15 

0. 

0. 

H+N02-N0+0H 

. 350E15 

0. 

1500. 

0+N02-N0+02 

.  100E14 

0. 

600. 

0+N2-N0+N 

.  184E15 

0. 

76250. 

N+O2-N0+0 

.  640E10 

1. 

6280. 

H+NO-N+OH 

. 222E15 

0. 

50500. 

HN0+HN0-N20+H20 

.  395E13 

0. 

5000. 

HN0+N0-N20+0H 

.  200E13 

0. 

26000. 

H+H02-H2+02 

.  251E14 

0. 

700. 

H+H02-OH+OH 

.251E15 

0. 

1900. 

H+H02-H20+0 

.  501E14 

0. 

1000. 

H02+0H-H20+02 

.  501E14 

0. 

1000. 

H02+0-02+0H 

. 501E14 

0. 

1000. 
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H02+H02-H202+02 

. 100E14 

0. 

1000. 

H202+0H-H20+H02 

. 100E14 

0. 

1800. 

H202+H-H20+OH 

. 316E15 

0. 

8940. 

H202+H-H02+H2 

. 170E13 

0. 

3750. 

CO+HO2-C02+OH 

. 575E14 

0. 

22930. 

CH20+H02-HC0+H202 

. 100E13 

0. 

8000. 

HC0+O2-C0+HO2 

. 331E13 

0. 

7000. 

HCN+OH-CN+H20 

. 440E13 

0. 

9000. 

HCN+O-NCO+H 

. 121E05 

2.64 

4980. 

HCN+ONH+CO 

. 517E04 

2.64 

4980. 

HCN+O-CN+OH 

.270E10 

1.58 

26600. 

CN+H2-HCN+H 

. S4SE12 

0.7 

4885. 

CN+O-CO+N 

.  180E14 

0. 

0. 

CN+02-NCO+0 

. 560E13 

0. 

0. 

CN+OH-NCO+H 

. 500E14 

0. 

0. 

NCO+H-NH+CO 

.  500E14 

0. 

0. 

NCO+O-NO+CO 

. 300E14 

0. 

0. 

NC0+N-N2+C0 

.  200E14 

0. 

0. 

NC0+N0-N20+C0 

.  190E14 

0. 

0. 

NC0+H2-HNC0+H 

.  858E13 

0. 

9000. 

HNC0+H-NH2+C0 

.  200E14 

0. 

3000. 

NH+H-N+H2 

.  300E14 

0. 

0. 

H+0H+M-H20+M 

.  750E24 

-2.6 

0. 

O+H+M-OH+M 

.  100E17 

0. 

0. 

H+02+M-H02+M 

»151E16 

0. 

-1000. 

H2CN+M-HCN+H+M 

. 100E17 

0. 

30000. 
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APPENDIX!  EXTENDED  MILLER-BOWMAN  MECHANISM  [11,30] 


Rate  constants  for  reactions  listed  in  form  k=AT'Bexp(-E/RT) 


Reaction 

A 

E 

CH4+OCH3+OH 

1.02E9 

1.5 

8604. 

CH4+02-CH3+H02 

7.90E13 

0.0 

56000. 

2CH3 <+M) -C2H6 (+M) 

LOW/1 . 7E38-6 . 051292 . / 
TROE/O. 6046927. 132./ 
H2/2/C0/2/C02/3/H20/5/ 

9.03E16 

-1.2 

654. 

CH3+H(+M)-CH4(+M) 

LOW/8. 00E26-3. 00.0/ 

SRI/0.45797.979./ 

H2/2/C0/2/C02/3/H20/5/ 

6.0E16 

-1.0 

0.0 

CH4+H-CH3+H2 

2.20E4 

3.000 

8750. 

CH4+0H-CH3+H20 

1.60E6 

2.1 

2460. 

CH4+H02-CH3+H202 

1.80E11 

0.0 

18700. 

CH3+H02-CH30+0H 

2.00E13 

0.0 

0.0 

CH3+02-CH30+0 

2.05E19 

-1.570 

29229. 

CH3+0-CH20+H 

8.00E13 

0.0 

0.0 

CH20H+H-CH3+0H 

1.00E14 

0.0 

0.0 

CH30+H-CH3+0H 

1.00E14 

0.0 

0.0 

CH3+0H-CH2+H20 

7.50E6 

2.0 

5000. 

CH3+H-CH2+H2 

9.00E13 

0.0 

15100. 

CH30+M-CH20+H+M 

1.00E14 

0.0 

25000. 

CH20H+M-CH20+H+M 

1.00E14 

0.0 

25000. 

CH30+H-CH20+H2 

2.00E13 

0.0 

0.0 

CH20H+H-CH20+H2 

2.00E13 

0.0 

0.0 

CH30+0H-CH20+H20 

1.00E13 

0.0 

0.0 

CH20H+0H-CH20+H20 

1.00E13 

0.0 

0.0 

CH30+0-CH20+OH 

1.00E13 

0.0 

0.0 

CH20H+O-CH20+OH 

1.00E13 

0.0 

0.0 

CH30+02-CH20+H02 

6.30E10 

0.0 

2600. 

CH20H+02-CH20+H02 

1.48E13 

0.0 

1500. 

CH2+H-CH+H2 

1.00E18 

-1.560 

0.0 

CH2+0H-CH+H20 

1.13E7 

2.0 

3000. 

CH2+OH-CH20+H 

2.50E13 

0.0 

0.0 

CH+02-HC0+O 

3.30E13 

0.0 

0.0 

CH+O-CO+H 

5.70E13 

0.0 

0.0 

CH+OH-HCO+H 

3.00E13 

0.0 

0.0 

CH+C02-HCO+CO 

3.40E12 

0.0 

690. 

CH+H-C+H2 

1.50E14 

0.0 

0.0 

CH+H20-CH20+H 

5.72E12 

0.0 

-380.0 

CH+CH20-CH2C0+H 

9.46E13 

0.0 

-515. 

CH+C2H2-C3H2+H 

1.00E14 

0.0 

0.0 

CH+CH2-C2H2+H 

4.00E13 

0.0 

0.0 

CH+CH3-C2H3+H 

3.00E13 

0.0 

0.0 

CH+CH4-C2H4+H 

6.00E13 

0.0 

0.0 

C+O2-C0+O 

2.00E13 

0.0 

0.0 

C+OH-CO+H 

5.00E13 

0.0 

0.0 

C+CH3-C2H2+H 

5.00E13 

0.0 

0.0 

C+CH2-C2H+H 

5.00E13 

0.0 

0.0 

CH2+C02-CH20+C0 

1.10E11 

0.0 

1000. 

CH2+0-C0+2H 

5.00E13 

0.0 

0.0 

CH2+0-C0+H2 

3.00E13 

0.0 

0.0 

CH2+02-C02+2H 

1.60E12 

0.0 

1000. 

CH2+O2-CH20+O 

5.00E13 

0.0 

9000. 

CH2+02-C02+H2 

6.90E11 

0.0 

500. 

CH2+02-C0+H20 

1.90E10 

0.0 

-1000. 

CH2+02-C0+0H+H 

8.60E10 

0.0 

-500. 

CH2+02-HC0+0H 

4.30E10 

0.0 

-500. 

CH20+0H-HC0+H20 

3.43E9 

1.180 

-447. 

CH20+H-HC0+H2 

2.19E8 

1.770 

3000. 

CH20+M-HC0+H+M 

3.31E16 

0.0 

81000. 

CH20+O-HC0+OH 

1.80E13 

0.0 

3080. 

HC0+0H-H20+C0 

1.00E14 

0.0 

0.0 

HCO+H-H+CO+M 

2.50E14 

0.0 

16802. 

CO/1 . 9/H2/1 . 9/CH4/2 . 8/C02/3 

.  0/H2O/5 . 0/ 

HCO+H-CO+H2 

1.19E13 

0.250 

0.0 

HCO+O-CO+OH 

3.00E13 

0.0 

0.0 

HCO+OC02+H 

3.00E13 

0.0 

0.0 

HC0+O2-HO2+C0 

3.30E13 

-0.400 

0.0 

C0+0+M-C02+M 

6.17E14 

0.0 

3000. 

C0+0H-C02+H 

1.51E7 

1.300 

-758. 

C0+02-C02+0 

1.60E13 

0.0 

41000. 

H02+C0-C02+0H 

5.80E13 

0.0 

22934. 

C2H6+CH3-C2H5+CH4 

5.50E-1 

4.000 

8300. 

C2H6+H-C2H5+H2 

5.40E2 

3.5 

5210. 

C2H6+0-C2H5+0H 

3.00E7 

2.000 

5115. 

C2H6+0H-C2H5+H20 

8.70E9 

1.05 

1810. 

C2H4+H-C2H3+H2 

1.10E14 

0.0 

8500. 

C2H4+0-CH3+HC0 

1.60E9 

1.200 

746. 

C2H4+0H-C2H3+H20 

2.02E13 

0.0 

5955. 

CH2+CH3-C2H4+H 

3.00E13 

0.0 

0.0 

H+C2H4 (+M) -C2H5 (+M) 

2.21E13 

0.0 

2066. 

LOW/6. 37E27-2. 8-54./ 

H2/2/C0/2/C02/3/H20/5/ 

C2H5+H-2CH3 

1.00E14 

0.0 

0.0 

C2H5+02-C2H4+H02 

8.43E11 

0.0 

3875. 

C2H2+0-CH2+C0 

1.02E7 

2.000 

1900. 

C2H2+0-HCC0+H 

1.02E7 

2.000 

1900. 

H2+C2H-C2H2+H 

4.09E5 

2.390 

864. 

H+C2H2 (+M) -C2H3 (+M) 

5.54E12 

0.0 

2410. 

LOW/2 . 67E27-3 . 52410 . / 

H2/2/C0/2/C02/3/H20/5/ 

C2H3+H-C2H2+H2 

4.00E13 

0.0 

0.0 

C2H3+0-CH2C0+H 

3.00E13 

0.0 

0.0 

C2H3+02-CH20+HC0 

4.00E12 

0.0 

-250. 

C2H3+0H-C2H2+H20 

5.00E12 

0.0 

0.0 

C2H3+CH2-C2H2+CH3 

3.00E13 

0.0 

0.0 

C2H3+C2H-2C2H2 

3.00E13 

0.0 

0.0 

C2H3+CH-CH2 +C2 H2 

5.00E13 

0.0 

0.0 

0H+C2H2-C2H+H20 

3.37E7 

2.000 

14000. 

0H+C2H2-HCC0H+H 

5.04E5 

2.300 

13500. 

OH+C2H2 “CH2CO+H 

2.18E-4 

4.500 

-1000. 

OH+C2H2-CH3+CO 

4.83E-4 

4.000 

-2000. 

HCCOH+H-CH2  CO+H 

1.00E13 

0.0 

0.0 

C2H2+0-C2H+0H 

3.16E15 

-0.6 

15000. 

CH2C0+0-C02+CH2 

1.75E12 

0.0 

1350. 

CH2CO+H-CH3+CO 

1.13E13 

0.0 

3428. 

CH2CO+H-HCCO+H2 

5.00E13 

0.0 

8000. 

CH2C0+0-HCC0+0H 

1.00E13 

0.0 

8000. 

CH2CO+OH-HCCO+H20 

7.50E12 

0.0 

2000. 

CH2CO  <+M) -CH2+C0 (+M) 

3.00E14 

0.0 

70980. 
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LOW/3. 60E150. 059270./ 
C2H+O2-2C0+H 

5.00E13 

0.0 

1500. 

C2H+C2H2-C4H2+H 

3.00E13 

0.0 

0.0 

H+HCCO-CH2 (S) +CO 

1.00E14 

0.0 

0.0 

0+HCC0-H+2C0 

1.00E14 

0.0 

0.0 

HCC0+02-2C0+0H 

1.60E12 

0.0 

854. 

CH+HCCO-C2H2+CO 

5.00E13 

0.0 

0.0 

2HCCO-C2H2+2CO 

1.00E13 

0.0 

0.0 

CH2 (S) +M-CH2+M 

1.00E13 

0.0 

0.0 

H/0.0/ 

CH2 (S)+CH4-2CH3 

4.00E13 

0.0 

0.0 

CH2 (S) +C2H6-CH3+C2H5 

1.20E14 

0.0 

0.0 

CH2 (S) +02— CO+OH+H 

3.00E13 

0.0 

0.0 

CH2 (S) +H2— CH3+H 

7.00E13 

0.0 

0.0 

CH2 (S) +H-CH2+H 

2.00E14 

0.0 

0.0 

C2  H+O— CH+CO 

5.00E13 

0.0 

0.0 

C2H+OH— HCCO+H 

2.00E13 

0.0 

0.0 

2CH2-C2H2+H2 

4.00E13 

0.0 

0.0 

CH2  +HCCO-C2  H3  +CO 

3.00E13 

0.0 

0.0 

CH2+C2H2-C3H3+H 

1.20E13 

0.0 

6600. 

C4H2+OH— C3H2+HCO 

6.66E12 

0.0 

-410. 

C3H2+02— HCO+HCCO 

1.00E13 

0.0 

0.0 

C3H3+02— CH2CO+HCO 

3.00E10 

0.0 

2868. 

C3H3+0— CH20+C2H 

2.00E13 

0.0 

0.0 

C3H3+OH— C3H2+H20 

2.00E13 

0.0 

0.0 

2C2H2— C4H3+H 

2.00E12 

0.0 

45900. 

C4H3+M— C4H2+H+M 

1.00E16 

0.0 

59700. 

CH2 (S) +C2H2-C3H3+H 

3.00E13 

0.0 

0.0 

C4H2+0— C3H2+CO 

1.20E12 

0.0 

0.0 

C2H2 +02— HCCO+OH 

2.00E8 

1.500 

30100. 

C  2  H  2 +M— C2 H+H+M 

4.20E16 

0.0 

107000. 

C2H4+M— C2H2+H2+M 

1.50E15 

0.0 

55800. 

C2H4+M— C2H3+H+M 

1.40E15 

0.0 

82360. 

H2+02-20H 

1.70E13 

0.0 

47780. 

OH+H2— H20+H 

1.17E9 

1.3 

3626. 

0+0H-02+H 

4.00E14 

-0.5 

0.0 

0+H2-0H+H 

5.06E4 

2.67 

6290. 

H+02+M-H02+M 

3.61E17 

-0.72 

0.0 

H20/18 . 6/C02/4 .2/H2/2 . 9/CO/2 
0H+H02— H20+02 

.  1/N2/1 .3/ 
•7.50E12 

0.0 

0.0 

H+H02-20H 

1.40E14 

0.0 

1073. 

0+H02-02+OH 

1.40E13 

0.0 

1073. 

20H-0+H20 

6.00E8 

1.3 

0.0 

2H+M-H2+M 

1.00E18 

-1.00 

0.0 

H2/0.0/H20/0.0/C02/0 .0/ 
2H+H2-2H2 

9.20E16 

-0.600 

0.0 

2H+H20-H2+H20 

6.00E19 

-1.25 

0.0 

2H+C02-H2+C02 

5.49E20 

-2.00 

0.0 

H+OH+M— H20+M 

1.60E22 

-2.00 

0.0 

H2O/5.0/ 

H+O+M-OH+M 

6.20E16 

1 

o 

o> 

o 

0.0 

H2O/5.0/ 

20+M-02+M 

1.89E13 

0.0 

-1788. 

H+H02-H2+02 

1.25E13 

0.0 

0.0 

2H02— H202+O2 

2.00E12 

0.0 

0.0 

H202+M-20H+M 

1.30E17 

0.0 

45500. 

H202+H-H02+H2 

1.60E12 

0.0 

3800. 

H202+0H-H20+H02 

1.00E13 

0.0 

1800. 
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CH+N2-HCN+N 

3.00E11 

0.0 

13600. 

CN+N-C+N2 

1.04E15 

-0.5 

0.0 

CH2+N2-HCN+NH 

1.00E13 

0.0 

74000. 

H2CN+N-N2+CH2 

2.00B13 

0.0 

0.0 

H2CN+M-HCN+H+M 

3.00E14 

0.0 

22000. 

C+NO-CN+O 

6.60E13 

0.0 

0.0 

CH+NO-HCN+O 

1.10E14 

0.0 

0.0 

CH2+NO-HCNO+H 

1.39E12 

0.0 

-1100. 

CH3+N0-HCN+H20 

1.00E11 

0.0 

15000. 

CH3+NO-H2CN+OH 

1.00E11 

0.0 

15000. 

HCCO+NO-HCNO+CO 

2.00E13 

0.0 

0.0 

CH2 (S) +NO-HCN+OH 

2.00E13 

0.0 

0.0 

HCNO+H-HCN+OH 

1.00E14 

0.0 

12000. 

CH2+N-HCN+H 

5.00E13 

0.0 

0.0 

CH+N-CN+H 

1.30E13 

0.0 

0.0 

C02+M-NO+CO 

1.90E11 

0.0 

3400. 

HCCO+N-HCN+CO 

5.00E13 

0.0 

0.0 

CH3+N-H2CN+H 

3.00E13 

0.0 

0.0 

C2H3+N-HCN+CH2 

2.00E13 

0.0 

0.0 

C3H3+N-HCN+C2H2 

1.00E13 

0.0 

0.0 

HCN+0H-CN+H20 

1.45E13 

0.0 

10929. 

OH+HCN-HOCN+H 

5.8SE4 

2.4 

12500. 

OH+HCN-HNCO+H 

1.98E-3 

4.0 

1000. 

OH+HCN-NH2+CO 

7.83E-4 

4.00 

4000. 

HOCN+H-HNCO+H 

1.00E13 

0.0 

0.0 

HCN+O-NCO+H 

1.38E4 

2.64 

4980. 

HCN+O-NH+CO 

3.45E3 

2.64 

4980. 

HCN+O-CN+OH 

2.70E9 

1.58 

26600. 

CN+H2-HCN+H 

2.95E5 

2.45 

2237. 

CN+O-CO+N 

1.80E13 

0.0 

0.0 

CN+O2-NC0+O 

5.60E12 

0.0 

0.0 

CN+OH-NCO+H 

6.00E13 

0.0 

0.0 

CN+HCN— C2N2+H 

2.00E13 

0.0 

0.0 

CN+N02-NC0+N0 

3.00E13 

0.0 

0.0 

CN+N20-NC0+N2 

1.00E13 

0.0 

0.0 

C2N2+0-NC0+CN 

4 . 57E12 

0.0 

8880. 

C2N2+OH-HOCN+CN 

1 . 86E11 

0.0 

2900. 

H02+N0-N02+0H 

2.11E12 

0.0 

-479. 

N02+H-N0+0H 

3.50E14 

0.0 

1500. 

N02+0-N0+02 

1.00E13 

0.0 

600. 

N02+M-N0+0+M 

1.1E16 

0.0 

66000. 

NCO+H-NH+CO 

5.00E13 

0.0 

0.0 

NCO+O-NO+CO 

2.00E13 

0.0 

0.0 

NCO+N-N2+CO 

2.00E13 

0.0 

0.0 

KCO+OH-NO+CO+H 

1.00E13 

0.0 

0.0 

NCO+M-N+CO+M 

3.10E16 

-0.5 

48000. 

NC0+N0-N20+C0 

1.00E13 

0.0 

-390. 

NCO+H2-HNCO+H 

8.58E12 

0.0 

9000. 

HNCO+H-NH2+CO 

2.00E13 

0.0 

3000. 

NH+02-HN0+0 

1.00E13 

0.0 

12000. 

NH+O2-N0+0H 

7.60E10 

0.0 

1530. 

NH+H0-N20+H 

2.40E1S 

-0.80 

0.0 

H20+0H-N2+H02 

2.00E12 

0.0 

10000. 

N20+H-N2+OH 

7.60E13 

0.0 

15200. 

N20+M-N2+0+M 

1.60E14 

0.0 

51600. 

N20+0-N2+02 

1.00E14 

0.0 

28200. 

N20+0-2N0 

1.00E14 

0.0 

28200. 

NH+OH-HNO+H 

2.00E13 

0.0 

0.0 
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NH+OH-N+H20  5 

NH+N-N2+H  3 

MH+H-N+H2  1 

NH2+0-HN0+H  6 

NH2+0-NH+OH  6 

NH2+0H-NH+H20  4 

NH2+H-NH+H2  6 

NH2+N0-NNH+0H  6 

NH2+N0-N2+H20  6 

NH3+OH-NH2+H20  2 

NH3+H-NH2+H2  6 

NH3+0-NH2+0H  2 

NNH-N2+H  1 

NNH+N0-N2+HN0  5 

NNH+H-N2+H2  1 

NNH+0H-N2+H20  5 

NNH+NH2-N2+NH3  5 

NNH+NH-N2+NH2  5 

NNH+ON20+H  1 

HNO+M-H+NO+M  1 

H2O/10 . 0/02/2 . 0/N2/2 . 0/H2/2 . 0/ 
HNO+OH-NO+H20  3 

HNO+H-H2+NO  5 

HNO+NH2-NH3+NO  2 

N+N0-N2+0  3 

N+02-N0+0  6 

N+OH-NO+H  3 

CH4+N02-CH3+H0N0  1 

CH3+N02-CH30+N0  1 

H02+N02-H0N0+02  4 

CH+N02-HC0+N0  1 

HONO+M-OH+NO+M  1 

N02+N02-N0+N0+02  2 

N02+N02-NO+N03  3 

HN03+M-N02+0H+M  1 

NO2+NO-N20+O2  1 

NO+M-N+O+M  4 

N0+N0-N2+O2  1 

NO+HNO-N20+OH  2 

N0+H02-HN0+02  2 

HC0+N02-C0+H0N0  1 

HC0+N02-C02+H+N0  1 

C0+N02-C02+N0  1 

H0N0+H-H2+N02  1 
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00E13 
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00E14 
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75E12 

0.0 

0.0 

00E6 

2.00 
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92E13 

0.0 

3650. 

40E15 

-1.25 

0.0 

20E15 

-1.25 

0.0 

04E6 

2.04 

566. 

36E5 

2.39 

10171. 

10E13 

0.0 

9000. 

00E4 

0.0 

0.0 

00E13 

0.0 

0.0 

00E14 

0.0 

0.0 

00E13 

0.0 

0.0 

00E13 

0.0 

0.0 

00E13 

0.0 

0.0 

00E14 

0.0 

0.0 

50E16 

0.0 

48680. 

60E13 

0.0 

0.0 

00E12 

0.0 

0.0 

00E13 

0.0 

1000. 

27E12 

0.30 

0.0 

40E9 

1.00 

6280. 

80E13 

0.0 

0.0 

20E14 

0.0 

30000. 

50E13 

0.0 

0.0 

64E11 

0.0 

-479. 

01E14 

0.0 

0.0 

8E17 

0.0 

44900. 

0E12 

0.0 

26825. 

9E11 

0.0 

23844. 

6E15 

0.0 

30600. 

0E12 

0.0 

60000. 

0E20 

-1.5 

150000. 

3E14 

0.0 

75630. 

0E12 

0.0 

26000. 

0E11 

0.0 

2000. 

1 5E13 

0.0 

-430. 

,  5E11 

0.0 
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26E14 

0.0 
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0E12 

0.0 
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